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Abstract 

That review article is focused on the tremendous progresses reahzed during the last 
fifteen years in the understanding of multifragmentation and its relationship to the 
liquid-gas phase diagram of nuclei and nuclear matter. The explosion of the whole 
nucleus, early predicted by Niels Bohr [1], is a very complex and rich subject which 
continues to fascinate nuclear physicists as well as theoreticians who extend the 
thermodynamics of phase transitions to finite systems. 
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1 Introduction 

The nuclear multifragmentation phenomenon was predicted and studied since 
the early 80 's. It is however only with the advent of powerful 47r detectors 
that real advances were made. Such arrays allow the detection of a large 
amount of the many fragments and light particles produced in nuclear colli- 
sions at intermediate and high energies. Indeed it now appears that further 
progresses are linked to the knowledge of many observables and the possi- 
bility to study correlations inside the multifragment events. This paper is a 
snapshot of what is known at the present time, making use of a large panel of 
the available data from heavy-ion collisions at intermediate energies, and from 
hadron-nucleus collisions. A comparison between fragmentation properties of 
quasi-projectiles formed in semi-peripheral collisions and of "fused" systems 
produced in the more violent ones will be particularly developed. The forma- 
tion of the fragmenting systems (or sources) in the course of nuclear collisions 
makes it necessary to model the collisions and several transport codes have 
been developed for two decades 0, 0], which account for many properties ex- 
perimentally observed. On the other hand statistical descriptions based on 
phase space occupation also account well for the static properties (partitions) 
of multifragmenting systems. The connection between both descriptions will 
be discussed. The Equation of State describing nuclear matter, similar to the 
van der Waals equation for classical fluids, foresees the existence of a liquid-gas 
type phase transition; multifragmentation was long-assimilated to this tran- 
sition ( or to a "liquid- fog" transition [sl]). Nuclear physicists are however 
dealing with finite systems - nuclei feeling nuclear and Coulomb forces - and 
not with infinite nuclear matter. Following the concepts of statistical physics, 
a new definition of phase transitions for such systems was recently proposed, 
showing that specific phase transition signatures could be expected. Different 
and coherent signals of phase transition have indeed been evidenced in a few 
cases. 

After a general survey of the multifragmentation phenomenon in section [21 
the necessity of sorting events and the different ways to proceed are pre- 
sented in section [31 The properties of the emitted fragments are detailed in 
section [H Section [51 introduces statistical and dynamical models whose re- 
sults are commonly compared to experimental data. The reconstruction of 
the multifragmenting systems connected with calorimetric and thermometric 
measurements (section ^ is followed by a study of the properties of these 
systems at the freeze-out stage in section [71 Finally the view of multifragmen- 
tation in terms of the phase transition of a finite system and the experimental 
signatures evidencing the transition are developed in sections M and M 
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2 A general survey 



One can come to the multifragmentation concept and its relation to a phase 
transition from two different starting points. Firstly by using kinetic models 
which show that nucleus-nucleus collisions at intermediate energies produce 
matter at subnormal density which breaks and where a phase transition is 
predicted to occur; it is a pure theoretical starting point. The other one con- 
sists in studying both experimentally and theoretically the evolution of decay 
mechanisms with increasing excitation energy. 



2.1 Nuclear matter at subnormal density and phase transition 



2.1.1 The nuclear liquid-gas phase transition 

Nuclear matter is an idealized macroscopic system with an equal number of 
neutrons and protons. It interacts via nuclear forces, and Coulomb forces are 
ignored due to its size. Its density p is spatially uniform. The nucleon-nucleon 
interaction is constituted by two components according to their radial in- 
terdistance : a very short-range repulsive part which takes into account the 
compressibility of the medium and a long-range attractive part. Changed by 
five orders of magnitude the nuclear interaction is similar to van der Waals 
forces acting in molecular medium. In a sense the phase transition in nuclear 
matter resembles the liquid-gas phase transition in classical fluids. However, 
as compared to classical fluids the main difference comes from the gas compo- 
sition: for nuclear matter the gas phase is predicted to be composed not only 
of single nucleons, neutrons and protons, but also of complex particles and 
fragments depending on temperature conditions 0, 0] ■ 

A set of isotherms for an equation of state (pressure versus density) corre- 
sponding to nuclear forces (Skyrme effective interaction and finite tempera- 
ture Hartree-Fock theory ^]) is shown in figure [H It exhibits the maximum- 
minimum structure typical of van der Waals equation of state. Depending on 
the effective interaction chosen and on the model j^, 0, [lol, 11|, the nuclear 



equation of state exhibits a critical point at pc ~-0.3-0.4po and ^16-18 MeV. 
Po and Vo refer to normal density /volume. The region below the dotted line 
in figure [1] corresponds to a domain of negative compressibility: at constant 
temperature an increase of density is associated to a decrease of pressure. 
Therefore in this region a single homogeneous phase is unstable and the sys- 
tem breaks into a liquid phase and a gas phase in equilibrium. It is the so 
called spinodal region, and spinodal instability corresponds to the breaking 
into the two phases. Such instability has been proposed, for a longtime, as 



a possible mechanism responsible for multifragmentation [12|, ll3|, [ij]. It will 



be discussed in section El The spinodal region constitutes the major part of 
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Fig. 1. Equation of state relating the pressure (left) or the temperature (right) and 
the density (normalised to critical values) in nuclear matter. The curves represent 
isotherms (left) and isobars (right). The dashed-dotted lines are the coexistence 
lines and the dotted lines the spinodal lines (from [ist]). 

the coexistence region (dashed-dotted line in figure [1]) which also contains two 
metastable regions: one at density below pc for the nucleation of drops and 
one above pc for the nucleation of bubbles (cavitation). 



2.1.2 From nuclear matter to hot nuclei 

Evidently the hot piece of nuclear matter produced in any nuclear collision 
has at more a few hundreds nucleons and so is not adequately described by 
the properties of infinite nuclear matter; surface and Coulomb effects can not 
be ignored. These effects have been evaluated and lead to a sizeable reduction 
of the critical temperature js], [q], [i6[. Finite size effects have been found to 
reduce the temperature by 2-6 MeV depending on the size of nuclei while 
the Coulomb force is responsible for a further reduction of 1-3 MeV. However 
large reductions due to small sizes are associated with small reductions from 
Coulomb. Consequently, in the range A=50-400 a total reduction of about 
7 MeV is calculated leading to a "critical" temperature of about 10 MeV for 
nuclei or hot pieces of nuclear matter produced in collisions between very heavy 
nuclei. The authors of reference ^ indicate that, due to some approximations, 
the derived values can be regarded as upper limits. Finally we can recall that, 
in infinite nuclear matter, the binding energy per particle is 16 MeV whereas 
it is about 8 MeV in a finite nucleus. Clearly these values well compare with 
the Tc values for infinite nuclear matter and nuclei just discussed. 
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2.2 Evolution of decay mechanisms: from evaporation to vaporization 



2. 2. 1 Evaporation 

The behaviour of nuclei at excitation energies around 1 MeV/nucleon has been 
extensively studied and rather well understood using statistical models [l7| . 
The theoretical treatment of particle emission involves the estimation of mi- 
crostate densities defined for equilibrium states. Such a treatment is justified 
only when there is enough time between successive emissions for the relax- 
ation of the emitting nucleus to a new equilibrium state. At such excitation 
energies the density stays very close to normal density of cold nuclear matter 
and the earliest evaporation model rests on the basic idea: an emitted parti- 
cle can be considered as originally situated somewhere on the surface of the 
emitting nucleus with a randomly directed velocity, it is why we use the term 
evaporation. Moreover particles are emitted sequentially and independently 
without any correlation. Note that in this excitation energy domain fission is 
also a deexcitation mechanism and can compete with evaporation. 



2.2.2 The multifragmentation domain 

At excitation energies comparable with the binding energy of nuclei a copious 
emission of particles and fragments is experimentally observed and the name 
"multifragmentation" was introduced. A first attempt to connect multifrag- 
mentation to a possible phase transition is to derive the excitation energy 
domain where this type of deexcitation is experimentally observed. Figure [2] 
displays, for different reactions, the evolution of the average fragment multi- 
plicity, normalized to the size of the multifragmenting system, as a function 
of the excitation energy per nucleon deposited into the system. All points fall 
on a single bell shape curve. The onset of multifragmentation takes place for 
excitation energies around 3 MeV/nucleon and the maximum for fragment 
production is found around 9 MeV/nucleon, i.e. close to the binding energy 
of nuclei. At higher excitation energy, due to the opening of the vaporiza- 
tion channel, the fragment production reduces. On the other hand, average 
time intervals between successive emissions have been estimated by analysing 
space-time correlations between fragments, taking advantage of proximity ef- 
fects induced by Coulomb repulsion. Figure [3] displays those average time in- 
tervals measured as a function of excitation energy deposited into the emitting 
system. A strong decrease of measured times with the increase of excitation 
energy is observed up to around 5 MeV/nucleon. Then a saturation appears 
around 50-100 fm/c which reflects the limit of sensitivity of the method. For 
such short times fragments can be considered as emitted quasi simultaneously 
and fragment emissions can not be treated independently. Note that, correla- 
tively, sequential statistical models fail in reproducing the observed emission 



6 



X 



I 



5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 
0.5 

0. 







' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' 

O QP Cl+Au 43 MeV/nucleon BeauUeu 
# QP Ge+Ti 35 MeV/nucleon Beaulieu 
▲ QP Au+Au ALADIN 

Xe+Sn central INDRA 
•k Ni+Au central INDRA 



o 

• A 

★ 



O 



1^ 

8 ° 
▲ 



o 



_s_l: 



o 



6 8 10 12 14 16 
Excitation energy (MeV/nucleon) 



Fig. 2. Average fragment multiplicity (normalized to the number of incident nucle- 
ons) as a function of the excitation energy per nucleon (data from 18, 19, i^, 2l|). 
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Fig. 3. Average fragment emission time as a function of the excitation energy per 
nucleon (data from 23, 24. [isl. jjl . 27 . 28, 29, [s^]- CC stands for central collisions, 
QP for quasi-projectile sources). 



2.2.3 Vaporization and identification of a gas phase 



At excitation energies around 10 MeV/nucleon, first experimental indications 
for the onset of vaporization (disinte grat ion into light particles with Z < 2) 
were reported in references (sil . Moreover a gas phase was iden- 

tified by comparison with a model by studying the deexcitation properties of 



vaporized quasi-projectiles produced in ^^Ar-|-^®Ni reactions 36|. Chemical 
composition and average kinetic energies of the different particles are well re- 
produced by a quantum statistical model (grandcanonical approach) describ- 
ing a real gas of fermions and bosons in thermal and chemical equilibrium. 
The evolution with excitation energy of the composition of vaporized quasi- 
projectiles is shown in figure HI Nucleon production increases with excitation 
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energy whereas emission of alpha particles, dominant at lower excitation ener- 
gies, strongly decreases. The regular behaviour observed is a strong indication 
that a rapid change of phase does not occur in the considered excitation en- 
ergy range. Note that an excluded volume correction due to finite particle 
size [s^] (van der Waals-like behaviour) was found decisive to obtain the ob- 
served agreement. In the model, the experimental range in excitation energy 




Fig. 4. Composition of vaporized quasi-pro- 
jectiles, formed in 95 MeV/nucleon 
36^j,_^58jYj collisions, as a function of their 
excitation energy per nucleon. Symbols are 
for data while the lines (dashed for He 
isotopes) are the results of the model. The 
temperature values used in the model are 
also given, (from [36] ) . 



Fig. 5. Average kinetic energies 
of particles emitted by vaporized 
quasi-projectiles at different ex- 
citation energies and formed in 
95 MeV/nucleon ^'^Ar+^^Ni colli- 
sions. Symbols are for data while full 
lines are the results of the model. The 
dashed lines refer to average kinetic 
energies of particles expected for an 
ideal gas (3T/2) at temperatures 
derived from the model, (from (s^). 



per nucleon of the source (9.5 to 27.5 MeV) was covered by varying the tem- 
perature from 10 to 25 MeV and the only free parameter, the free volume, was 
fixed at 3Vo in order to reproduce the experimental ratio between the proton 
and alpha yields at 18.5 MeV/nucleon excitation energy. The average kinetic 
energies of the different charged particles are also rather well reproduced over 
the whole excitation energy range (fig. [5]) but the model fails to accurately fol- 
low the dependence on the different species especially for alphas. The dashed 
lines in fig. [5] indicate the average kinetic energies, 3T/2, expected for an ideal 
gas, which appear as a rather good approximation; this is due to the low den- 
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sity of the emitting source. We are in presence of a quantum weakly-interacting 
gas. 



3 Experimental event sorting 

Multifragmentation is by essence associated to the emission of several frag- 
ments. Any study of the phenomenon requires a coincident and efficient detec- 
tion of these fragments and of the associated particles (Z<2). This is why, in 
the recent years, multifragmentation studies were performed with An detectors. 
These do not however guarantee a full detection in all reaction configurations. 
The combined geometrical efficiency (due to dead zones between detection 
cells and the necessary free space close to and 180° to allow beam propaga- 
tion) and detection/identification thresholds of charged particle (CP) arrays 
for instance, strongly disfavour any correct detection of the most peripheral 
collisions for energies above ~20 MeV/nucleon; indeed the projectile residues, 
emitted very close to 0°, escape detection, while the low energy of the target 
residues makes them unidentifiable, except if their time of fiight is measured. 
Moreover for heavy systems a very large number of neutrons is emitted before 
the release of any light charged particle as shown in fig. [61 Coupling a Air CP 
array to a neutron ball partly remedies this last drawback. 
For an efficient study of multifragmentation, the notion of complete event was 
proposed: a high enough part of the reaction products should be correctly de- 
tected and identified. For central collisions, a condition on the total detected 
charge (more than 70-90% of the system charge) is set. For quasi-projectile, 
one requires the detection, in the forward centre of mass hemisphere, of a large 
fraction of the projectile charge, or momentum. Owing to the non-detection 
of neutrons, a pseudo momentum, ZV, is sometimes used, replacing the mass 
of each product by its charge. In any case, the representative character of the 
selected events should be verified. 

Studies on multifragmentation should apply to homogeneous samples of events, 
which requires an appropriate sorting; this is mandatory for thermodynam- 
ical studies for instance. In hadron-nucleus collisions all events have similar 
topological properties independently of the impact parameter, as a single hot 
nucleus is created after a more or less abundant preequilibrium emission. Con- 
versely, in heavy-ion collisions, the outgoing channel is different depending on 
the masses and asymmetry of the incident partners, the incident energy and 
the impact parameter. At intermediate energies residual interactions (NN col- 
lisions) strongly compete with mean field effects; the number of NN collisions 
largely fluctuates, leading to different flnal reaction channels for the same 
initial conditions. The weakening of the mean fleld hinders, on average, full 
stopping above about 30 MeV/nucleon; the large fluctuations mentioned above 
allow however the observation of "fusion" (one outgoing heavy fragment) at 
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Fig. 6. Neutron vs light charged particle multiplicities for Xe+Bi reactions at var- 
ious incident energies. Larger multiplicities correspond to more violent collisions. 



From 39|] 



higher energies, although with small cross sections [38| . Most of the collisions 
end-up in two remnants coming from the projectile and the target - accompa- 
nied by some evaporated particles -, and some fragments and particles with 
velocities intermediate between those of the remnants: these are called mid- 
velocity products. They may have several origins, e.g. direct preequilibrium 
emission from the overlap region between the incident partners, or a neck of 
matter between them which may finally separate from QP or QT, or both. 
Whatever the type of reaction is, a fraction of the incident translational energy 
is transformed into "excitation energy" , E*, which may be shared into thermal 
energy (heat) and collective energies (rotational, expansion. . . ). While several 
experimental methods give access to E*, the knowledge of the fraction allot- 
ted to thermal or collective energies relies on models. The sorting of events, 
is generally done through global variables, which condensate the information 
measured on each event. Two philosophies guide the methods used for event 
sorting: the impact parameter dependence, and the event topology. 



3.1 Impact parameter selectors 



Sorting events with respect to an "experimental" impact parameter, bexp is 
appropriate for studying global properties of collisions versus their violence. 
These methods are also useful for comparisons with transport codes, pro- 
vided their outputs are filtered and sorted accordingly. Some observables are 
strongly connected to the impact parameter and are thus commonly used as 
sorting variables. One of the most popular impact parameter selector (IPS) 



is the charged product multiplicity [40|, sometimes reduced to the - barely 
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Fig. 7. Total transverse energy distributions for 
various impact parameters, calculated in the 
DYWAN simulation, for the 95 MeV/nucleon 
Ar+Ni system, (private communication from F. 
Sebille). 



smaller - light charged particle (Icp, Z=l,2) multiplicity. The fragment (Z>3) 
multiplicity is in no way a good selector, because its small value, less than 10, 
induces too large fluctuations. IPS based on coupled neutron and Icp multi- 
plicities are used with neutron balls 4l|, |42]. Other IPS are the total charge 
bound in charged products, excluding hydrogen isotopes, Zhound = > 2), 



Fig. 8. Relation between the 
dissipated energy, Ediss (see 
text), and the impact parameter 
for the Ni+Au reactions at 
52 MeV/nucleon modelled with a 
stochastic mean field simulation. 
From 4£ 



introduced by the ALADIN group |43 
forward emitted in the centre of mass 



or the value of the largest charge 



4J], Z^^^. This last variable requires 



arrays able to identify products in a very large range of atomic number and 
angles. Implicit in the previous IPS is the notion of dissipation, namely the 
part of the initial translational energy transformed into excitation energy of 
the system. Several IPS reflect more directly this connection, and first of all 
the transverse (perpendicular to the beam) component of the energy, either 
of all products [461], or of Icp only j47|. The latter choice is due to the gener- 
ally better efficiency of the arrays to Icp than to fragments. At intermediate 
energies, a narrow zone of b^xp generally covers a broad range of true impact 
parameters, as shown in fig. [71 obtained with the DYWAN code 48j. for the 
Ar+Ni system at 95 MeV/nucleon. If Z^^ is identified, its velocity, Vzav^, can 
be used to quantify dissipation [i^. A related variable is the energy dissipa- 
tion calculated as if the reaction was purely binary, without mass transfer [5cJ]: 
Ediss = Ecm — 1/2 /X VzAv X {Ap + At) /At- Fig. [H] shows for instance the good 
correlation between the impact parameter and Ediss, calculated in stochastic 
mean field simulations for the Ni-|-Au reactions at 52 MeV/nucleon, despite 
the presence of mid-rapidity emission in the simulation; in this figure stars and 



squares correspond to simulations with different isospin terms in the EOS [51 



Finally correlations between the studied variables and the sorting variables 
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Fig. 9. Means, variances and normalized variances of the fragment (Z>3, left) or 
IMF (3< Z < 20, right) multiplicities versus the beam energy for central collisions 
between ^"^^Xe and ""*Sn (INDRA data) and ^^^Xe and ^^^La (from using 
different IPS: multiplicity of charged products, Nc, of Icp, M/cp, of hydrogen, Nh, 
total detected charge, Zdet, charge bound in Icp, Zicp, transverse energy total, Etr, 
or of Icp, Ei2- 



must be avoided as much as possible. A detailed study on this subject was 



published in [52], where the average and the variance of the fragment mul- 
tiplicity distributions obtained when selecting central collisions with differ- 
ent IPS (upper 10% of the IPS distribution) were examined. The authors 
found that Mjmf was positively auto-correlated (small variance and normal- 
ized variance) with the total detected charge and negatively auto-correlated 
(small variance and large normalized variance) with variables related to Icp or 
proton multiplicity. They conclude that the more neutral selectors were the 
charged product multiplicity and transverse energy (fig. [9]). These conclusions 
on the best IPS are detector dependent: the same study, on a similar system 
but using the INDRA array, shows that the negative correlation with Icp re- 
lated variables weakly persists. The total detected charge (complete events) 
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is, for this system, an IPS as neutral as the charged product multiphcity and 
total transverse energy. 



3.2 Topology selectors 



IPS classify events with respect to the violence of the collisions, related to 
the energy transformed into thermal energy, irrespective of the kind of reac- 
tion which happened. Indeed, at intermediate energies (10-100 MeV/nucleon), 
there are large fluctuations in the exit channel associated to a given impact pa- 



rameter, as shown in stochastic transport models [53| . For this reason, and par- 
ticularly when it appeared desirable to work on a rather well defined "source" , 
sorting variables related to the topology of the events, namely to their shape 
in velocity space, were proposed. Most of the following variables are borrowed 
from high energy physics. 

The flow angle, 9 flow, is the angle between the beam axis and the main axis 
of the kinetic energy-flow tensor [54], Qij = J2n=i PnPn/'^''^n, with m„ the 
mass and the i^^ Cartesian component of the momentum of particle n. The 
eigenvectors of the tensor define an ellipsoid, which corresponds to a rotated 
reference frame with respect to the centre of mass. The flow angle is connected 



with the impact parameter [55[ . By analogy with the evolution of low energy 
deeply inelastic collisions, it is expected that this angle increases with the 
violence of the collision while the exit channel is dominantly binary; around 
the Fermi energy the rotation angle of the system remains small so the bulk 
of the cross section is associated to small flow angles. If "fusion" occurs no 
privileged direction is expected, so this type of events is better isolated for 
large flow angles. This selection was largely used in the INDRA collaboration 
for the study of compact single source events. 

Among shape global variables are the isotropy ratio, Riso, and the energy 



ratio 
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where 6'„ is the emission angle of particle n in the centre of mass (or in the el- 
hpsoid) frame, and £"„ its kinetic energy. Both Riso and E^at quantify the part 
of the momentum or energy transferred from the beam to the transverse direc- 
tion; they are thus equal to zero for peripheral collisions and take large values 
for the more central collisions, associated to more spherical shapes. Large flow 
angles, and Riso or E^at values close to 1 equally select events with compact 
shapes. A very powerful variable, which would deserve a larger utilisation, is 
the charge density, Pzik) = E Zi(k) / Y Zj, where k is the centre of mass ve- 
locity of particle i (with charge Zi) [57|. With this variable it is possible to 
isolate rare events which are either binary - without mid-rapidity emission - or 
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monosource, from the majority of events which exhibit a mid-rapidity compo- 
nent between quasi-projectile and quasi-target remnants. All these selections 
implicitly assume the isotropy - or at least the forward-backward symmetry - 
of any source emission. 



3.3 Multivariate analysis techniques 



Event selections with both IPS and topology selectors are made via a sharp 
cut-off in the distribution of one or two variables. While it has the advan- 
tage to be simple and transparent, the physical implications might be more 
difficult to evaluate due to the large fluctuations of most of the global vari- 
ables. A generalisation to more variables is given by the multivariate analysis 
techniques 58|]. Principal Component Analyses allow reduction of the dimen- 
sionality of the information, and are based on a set of global variables. The 
multidimensional space is rearranged on axes which are linear combinations of 
the initial global variables, and carry a more or less important part of the in- 
formation. Projecting the events on the plane defined by the two axes bearing 
maximum information allows to separate classes of events with a rather close 

aim 



60, 61, 62 



topology [2l|, |59(]. Discriminant Analyses, neuronal networks 
at discriminating types of events. They need a learning phase before being 
applied to a p hysical problem, which is often performed with the help of a 



model [63|, l6J] , but can as well directly utilize a sample of real events [65|, [66 
Indeed if the model used does not correctly consider and weigh all the possi- 
ble reaction channels, the final results may be biased. The advantage of these 
methods is that they generally provide samples of events with a higher statis- 
tics than the sharp-cut ones. The properties of the samples must however be 
studied in detail in order to be sure to have selected an homogeneous ensemble 
of events. 



3.4 Cross sections 



The knowledge of absolute cross sections is of great interest, whether one 
sorts events with an IPS or isolates one emission source. In the first case, 
for a comparison of experimental results with those of a transport code, the 
knowledge of the correct zone of impact parameter is mandatory. Indeed an 
experimental impact parameter, hexp-, is estimated with the help of a variable 
$ assumed to vary monotonously with 6, using a geometric prescription (67|: 
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In the above formula, $ decreases for increasing values of b, from a value ^max 
for 6=0 to for 6 = bmax- Obviously bexp can be compared to the true impact 
parameter only if bmax is associated to the reaction cross section, or if the total 
measured absolute cross section is known. In the case of the study of a source 
properties, the interest is to find out whether the observed phenomenon is 
dominant or rare. 



3. 5 Conclusion 



It was shown in this section that, when dealing with the results of power- 
ful An arrays, it is necessary to sort the data. Several methods are used by 
different groups, based on global variables which condensate the information. 
The adequate sorting differs depending on the type of analysis which is fore- 
seen, either an evolution with the violence of the collision and a comparison 
with dynamical simulations, or a thermodynamical study of the properties 
of a source. As an example, for central collisions, it is known that topology 
selectors and IPS do not select the same events, due to fiuctuations in the 
reaction mechanism associated to a given impact parameter; while the former 
isolate "compact shape events". IPS favour events with elongated shapes in 
velocity space. It must also be stressed that any sorting is detector dependent, 
even if the "same" global variable is adopted; this must be kept in mind when 
comparing data obtained with different experimental set-up. 



4 Fragment properties 



In this section several fragment properties will be discussed, either static like 
multiphcities or charge distributions, or dynamical. The "references" are cho- 
sen as raw as possible, namely the total system or the projectile masses, and 
the available centre of mass energy, to facilitate comparisons with dynamical 
simulations. The main trends displayed in this section are however represen- 
tative of those found for the properties of well isolated emitting sources versus 
the excitation energy, as it will appear in the next sections. 
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Fig. 10. Average fragment (Z>3) mul- 
tiplicities as a function of the centre of 
mass energy for central collisions and sys- 
tems of various sizes, measured with IN- 
DRA. Statistical error bars are smaller 
than the size of the symbols. Systematic 
errors are reported for the GSI data. 



Fig. 11. Maximum average fragment mul- 
tiplicities (bottom panel) and the cen- 
tre of mass energy of its occurrence (top 
panel) for systems of various sizes. Cen- 
tral collisions (stars) measured with IN- 
DRA and the MSU 4tt array; quasi-pro- 
jectiles (circles) measured by EOS, AL- 
ADIN and Laval [l^. The lines are lin- 
ear fits of central (full) and QP (dashed) 
data. 



4-1 Multiplicities 



The fragment multiplicities in multifragmentation were abundantly studied in 
the last years. The most prominent observation is a rise and fall of the average 
fragment multiplicity, observed since the early 90's for Au quasi-projectiles [68| 
and for Au+Au central collisions [Slj: below 100 MeV/nucleon the fragment 
multiplicity increases when going from peripheral to central collisions, while 
above that beam energy the maximum fragment multiplicity is observed in 
semi-central collisions. The absolute maximum value is however reached in 
central collisions at rather moderate energy. A first study of the evolution 
of the IMF (3<Z<20) multiplicity versus the bombardin g en ergy and the 



size of the system was performed for central collisions in [69|. In the same 
line, we show in fig. [TD] similar results obtained with the INDRA array op- 
erating at GANIL and GSI. Central collisions were selected, as in 6^, by a 



cut of the top 10% of the total transverse energy distribution; we chose here 
to present fragment multiplicities, Mf, fragments being defined as all fully 
identified products with a charge Z>3. The identification thresholds are re- 
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spectively 0.8, 1.1, 1.3, 1.7 MeV/nucleon for Z=3, 10, 20, 50. Lowering all 
thresholds below 0.8 MeV/nucleon for all species would increase the average 
multiplicities by less than 3% for symmetric systems and up to a maximum of 
10% for Ar+Ni at 32 MeV/nucleon. Larger multiplicities are found when the 
mass of the system increases. Systematic errors, essentially due to efficiency, 
can be estimated by requiring additionally a high completeness of the events; 
they lie between 0.5 and 1.5 units, increasing with the incident energy, and 
are larger for the GSI data due to the dysfunction of some INDRA modules. 
With increasing bombarding energy, the rise and fall of the average fragment 
multiplicity is visible in all cases (but perhaps for the lighter system Ar+KCl, 
where there is no rising part in the measured energy range); Mf rises rapidly 
at low energies up to its maximum value, M™"^^, but decreases smoothly on 
the high energy side; indeed sizeable fragment multiplicities persist in central 



collisions at much higher energies [70|, l7l|]. The cm. energy at which the frag 



ment multiplicity is the largest also increases with the size of the system. In 
order to quantify these effects, the data of fig. [10] were fitted with a third de- 
gree polynomial; for comparison, the data of |69] were fitted in the same way. 
Figure [11] shows, as a function of the total mass of the systems, the maximum 
fragment multiplicity, MJ*"^, and the energy at which this maximum occurs. 
The maximum multiplicities measured with the MSU 47r array and with IN- 
DRA exhibit a remarkably coherent behaviour, despite the differences in the 
detection devices: M^""^ has a value around 2 for masses 90-100, and grows 
further on proportionally to the total mass of the system (2.7 units of multi- 
plicity for 100 incident nucleons). The consistency of the results indicates that 
]\^max -g sensitive to the highest charge limit included in Mf - indeed this 
limit only affects the tail of the largest fragment distribution at low energies; 
conversely it strongly depends on the low charge threshold chosen, as shown 
in 



72l |. No saturation of MJ^"'^ at the highest masses, as quoted in [69], due 



to the Kr-|-Au data of [tII], is visible. The energy where MJ'"'^ is observed 
first rapidly increases with the system mass, then tends to level off beyond 
mass 150 for the INDRA data. Note that E'^^ is systematically higher for the 
MSU data, although its evolution with Agys is parallel to that of the INDRA 
data. is also sensitive to the low charge limit included in Mj, decreasing 

when this limit is raised. The increase of the degree of fragmentation with the 
available energy is expected if thermal or radial flow energy put in the system 
increases. When it becomes very high, more hydrogen and helium isotopes are 
formed at break-up, and during the deexcitation stage, leading to a fall of the 
fragment multiplicity. Fragments being smaller for lighter systems can be a 
flrst reason of the early drop of fragment multiplicity for such systems. 
Finally it is also instructive to compare the fragmentation of quasi-projectiles, 
as those studied by the Laval group, the ALADIN and EOS collaborations, 
with that of medium mass to heavy systems formed in central collisions, to 
get a flrst hint on the similarities or differences between the fragmentation 
mechanisms. The maximum QP fragment multiplicity (circles in flg. [TTj) rises 
linearly with the projectile mass; the values remain however smaller than those 
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Fig. 12. Normalized average fragment 
multiplicities versus normalized Zhound 
(Z>2) for central collisions measured 
with INDRA and for Au and U quasi-pro- 
jectiles detected with ALADIN (adapted 
from figs. 6 and 10 of ref. [19]) The hor- 
izontal error bars for the INDRA data 
represent the RMS of the Zhound distri- 
butions. 
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Fig. 13. Charge asymmetry of fragments 
as a function of the excitation energy of 
the emitting source, for Au quasi-pro- 
jectiles (open circles), monosources from 
central Gd+U collisions (triangles) and 
Xe+Sn from 25 to 50 MeV/nucleon 
(other symbols). From [t^I 



reached in central collisions. The values of MJ'"'^ 



for q uasi-pro i ectiles are in- 
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dependent of both the target mass and the incident energy 
indication that multifragmentation is mainly driven by the thermal energy 
deposited in the system, as noted in |75|]. The difference between the maxi- 
mum multiplicities for QP's and central collisions might then come from the 
larger expansion energy found in the latter (see next subsection). Figure 9 of 
ref. [l^ showed a scaling of the average fragment multiplicity versus Zbound 
when both quantities were normalized to the projectile charge (Xe, Au or U). 
This scaled result is summarized in fig. [T21 Zbound was in that case taken as 
an IPS. For central collision samples measured with INDRA and used in this 
section, the average multiplicity and Zbound (and its RMS) were normalized 
to the total charge of the system: a scaling is also observed for the Xe+Sn 
and Au+Au systems, but fails for Ni+Ni (Ztot=56), whereas it persists for Xe 
quasi- projectile (Zp=54) [1^. The fragment multiplicity- Zbo„„(i scaling is thus 
a property of multifragmentation and not just a geometrical property when 
Zbound is taken as an IPS. As already observed, for heavy systems, the scaled 
maximum multiplicity is higher in central collisions and occurs at higher scaled 
Zbound, meaning that the multifragmentation partitions are different in both 
types of collisions. This is confirmed in fig. [T31 where the average charge asym- 
metry of fragments, Az = az/{{Z) x ^/Mj — 1) is plotted versus the excitation 
energy of QP sources or monosources [7J| (see sect. [6] for the calculation of 
E*): the asymmetry is smaller for monosources formed in central collisions. 
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Xe+Sn, collisions centrales 




Fig. 14. Charge distributions measured 
for central Xe on Sn collisions between 
25 and 100 MeV/nucleon incident en- 
ergies. From [t?! 



Fig. 15. Average centre of mass ki- 
netic energy of fragments as a func- 
tion of their charge for central Xe 
on Sn collisions (selected with the 
Icp transverse energy) between 25 and 
100 MeV/nucleon incident energies. 
From [TTj] 



which means that the system is more fragmented. 



4-2 -^feownd ^''^d, charge distributions in central collisions 



While the fragment multiplicity exhibits a rise and fall with the incident en- 
ergy, the total charged product multiplicity keeps increasing, indicating that 
the formed fragments become smaller. Indeed the total charge bound in frag- 
mentSj__z£,^„^ = > 3), monotonously diminishes with increasing en- 

ergy 72|; for heavy systems (Xe+Sn, Au+Au) the Z^^^^^ distributions shift 
to smaller values and become narrower when raising the incident energy from 
30 to 150 MeV/nucleon. For the Ni+Ni system it appears that the high Z 
tail of these distributions barely evolves above 80 MeV/nucleon. The distri- 
butions are no longer symmetric, but raise for low Zl^^^^^, indicating that the 
system goes toward vaporization. The charge distributions show the same evo- 
lution, extending over a narrower charge range at higher energies, as shown 
in fig. [H] for the Xe+Sn system. The shape of the distribution comprises a 
broad plateau at moderate energies, and evolves towards an exponential de- 
crease above 100 MeV/nucleon (see also fz^l). The yield of fragments with 
charges Z=3-10 is almost independent of the incident energy, between 39 and 
150 MeV/nucleon, for this system. An interesting property of the charge dis- 
tributions is their independence, when scaled by the fragment multiplicity, 
with respect to the total charge of the (heavy) system (104 and 156 in |76|), 
provided that the available energy is similar. 
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4-3 Energy distributions - expansion 



Energy spectra in the cm. frame of fragments emitted by compact multifrag- 
menting sources were shown in [78] for the two heavy systems 32 MeV/nucleon 
Xe+Sn and 36 MeV/nucleon Gd+U. They have asymmetric shapes for the 
hghter fragments, and tend to become more symmetric for Z>15. These spec- 
tra bear information about the Coulomb repulsion, the radial expansion and 



the temperature of the source from which they originate. Indeed in [79| the 
Coulomb barrier deduced from the maximum of some fragment energy spectra 
was used to infer the density of the emitting source; a warning was made on 
this method by S^l, where it was shown that the same shift of the maximum 
of the spectrum may also be due to the total Coulomb energy which increases 
faster than the fragment multiplicity. In any case fragment energy spectra 



are precious probes to test dynamical multifragmentation models [78|, |8l| . For 



compact source selection, fragment emission is roughly isotropic in the cen- 
tre of mass. This is not the case for selections based on transverse energy. 



ACP [82| or neural networks [83|. In this case either the sources are shown 
to be elongated along the beam axis in coordinate space, and the fragment 
energy is larger in this direction, or the selection includes collisions with a 
binary character in the exit channel (see section [3]). 

The average kinetic energy of fragments as a function of their charge is thus 
dependent on the event selection and in some cases on the emission angle, 
except that of the lighter fragments which is more stable against experimental 
selections. For compact sources formed in central collisions between 25 and 
50 MeV/nucleon Xe and Sn, the average cm. kinetic energy is rising with the 
fragment charge, then saturates and even decreases for charges Z>20-25. The 
decrease is essentially due to the properties of the largest fragments of the 
partitions. This evolution, also observed in j8l|, pleads for the kinetic energy 
essentially originating from Coulomb (and possibly radial expansion) effects 
rather than from thermal motion; in this latter case indeed, the average energy 
would be constant whatever the fragment charge or mass. Conversely with a 
transverse energy selection, the kinetic energy continuously increases with Z 
at 50 MeV/nucleon. More information was derived in 78| from the average 
kinetic of fragments as a function of their charge, for events sorted according 
to the fragment multiplicity and the rank of the fragment in the event. For 
each of the two above-mentioned systems, the experimental pattern was found 
independent of the fragment multiplicity, indicating that all emitting sources 
have similar charges. The largest fragment displays a specific behaviour which 
will be discussed later. The kinetic energy, for a given charge, increases with 
the collision energy, (see figH^and fig.39 in Mi)- If the mass and the volume 
of the source (i.e. Coulomb effects) remain constant, then the extra energy 
has to come from an expansion which increases with the incident energy. 
For a long time dynamical models have predicted, for central collisions be- 
tween heavy nuclei, the occurrence of a compression phase followed by an 
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Fig. 16. Systematics of radial expansion energy per nucleon as a function of the 
available energy for central collisions (see text). Data are taken from: Ar+Sc 84i |. 
Zn+Au |8l|, Kr+Au, Au+Au MSU [sj, Au+Au MULTICS H, Pb+Au Nau- 
tilus fs?] and INDRA (see text). 



expansion of the system [88|, l89|; the total collective energy (thermal plus 
compressional effects) can be derived at different times during the collision 
and at different distances from the centre of the expanding system. This last 
point gives some information about the evolution of radial velocity flow as a 
function of the radial distance r. As for the "Big Bang" a self similar expansion 
(collective velocity proportional to r) is observed up to (~ 80-100 fm/c after 
the beginning of the collisions, when the system is still homogeneous (pre- 
fragments begin to appear at that time). With such a prescription the radial 
expansion velocity is extracted from experiments using: Vrarf(r) = {r/R)vo, 
where R is the rms of pre-fragment distances to centre and Vq the flow ve- 
locity. The radial expansion energy E{R) is related to the average expansion 
energy through: {Erad) = 3E{R)/5; vq and E{R) are used as reference values 
for comparisons. Note that a Coulomb repulsion when fragments are emit- 
ted from a spherical volume of uniform density also gives a collective velocity 
field proportional to r. Consequently the two effects can not be distinguished 
and radial expansions are extracted from experiments with the help of models 
as SMM or event generators as SIMON, assuming given volumes at freeze- 
out [qoI, 91]. When dealing with deformed sources, a different flow proflle was 



found, Vrad(r) = (r/i?)"t'o, with a varying between 1 and 2 [92|. Expansion 
energy is not a generic feature of statistical multifragmentation models, the 
implicit assumption is that partitions result only from the thermal part of the 
energy and are not influenced by the flow component. This hypothesis was 

even 



shown to be correct in the framework of the lattice gas model [9^ 
when the flow energy amounts to 50% of the available energy. 



94 



A systematics of radial expansion energies derived for central collisions and 
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various systems in the Fermi energy domain is reported in figure [T6[In most 
cases, it was derived from comparison witli a statistical model. In (8l| it was 
obtained by subtracting the kinetic energies obtained at low incident energy, 
as representative of the Coulomb component, from those measured at higher 
energy; for compact sources, the statistical model implies a spherical volume, 
while for other selections the authors used SMM with a deformed envelope, 
or consider only the energy at 90° in the cm. For the INDRA data quoted in 
fig- [m several determinations were made for the same systems. Average values 
are reported in the figure 2^, 65, 82, 91, 92, 95]. Despite the differences in the 



methods of determination of the radial energy, a general trend emerges, namely 
the onset of expansion energy around ~ 5 MeV/ nucleon followed by an increase 
with the cm. available energy. Values remain small (<1 MeV/nucleon) up 
to available energies around 8-9 MeV/nucleon. No clear size effect can be 
inferred from the available data, because of large error bars, and the difficulty 
to disentangle Coulomb and expansion energies. There may be an indication 
that there is less expansion when the entrance channel is very asymmetric, 
supported by the Ni+Au and Zn+Au data, but not by those on Kr+Au. 
Note also that radial expansion energies deduced from the BOB dynamical 
simulations for the Xe+Sn system at 32 MeV/nucleon and for the Gd+U 
system at 36 MeV/nucleon, both studied with INDRA, are in good agreement 
with this systematics (9]| . 

It is much more hazardous to derive fragment kinetic energies, and expansion 
energies, from quasi-projectiles, as in this case the reference frame is built 
from the fragment kinetic properties. The indications found in the literature 
are however that there is no expansion for QPs as long as their excitation 
energy, e*, does not reach ~5 MeV/nucleon; for light Ar QPs, the expansion 
energy is reported to increase from 0.2 to 0.7-1.1 MeV/nucleon for e* varying 
between 6 and 8 MeV/nucleon j96|, while for heavy Au-like nuclei an expansion 
of 0.8 MeV/nucleon was found at e* = 6 MeV/nucleon [97]. A similar small 
variation versus the excitation energy of the source was observed in hadron 
induced reactions on Au nuclei, with an expansion increasing from 0.2 to 
0.6 MeV/nucleon when e* varies between 5.5 and 8 MeV/nucleon [q^. In 
these two cases, the extra energy can be attributed to thermal pressure only, 
while for central collisions between heavy ions a compression phase occurs. The 
drawback on the reference frame for QP's can be overcome by using relative 
velocities between fragments, in order to compare values from the different 



kinds of multifragmenting sources [74 



4-4 The specific role of the largest fragment 



Several studies point out the specific properties of the largest fragment in each 
partition, be they static or dynamic. In the MMMC statistical model it was 
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Fig. 17. Average charge of the largest 
fragment versus the available energy 
for transverse energy selected central 
collisions measured with INDRA. 
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Fig. 18. Average charge of the largest frag- 
ment versus the total system mass for 
three available energies and central colli- 
sions. 



shown that the size of the largest fragment strongly decreases with increasing 
thermal excitation energy of the source but is independent of the source total 
charge (when varied between 56 and 90); it may thus be used to estimate 
this energy [99|]. The largest fragment size also decreases when the volume 
occupied by the source grows larger, as indicated in [lOQ], while those of the 
fragments of higher rank are roughly constant. It must be noted that such a 
variation of the ranked fragments was observed when plotting the fragment 
size as a function of the fragment multiplicity [78t|. 

From an experimental point of view, the independence of the largest frag- 
ment charge with respect to the system size follows from that of the charge 
distribution [t^ and was indeed observed for the 32 MeV/nucleon Xe+Sn 
and 36 MeV/nucleon Gd+U compact single sources isolated among central 
collisions. Such a study can be extended to the central collisions studied with 
INDRA already used in the previous subsections (transverse energy selection). 
The average charge of the largest fragment, for a given system, first strongly 
decreases with increasing available energy and then tends to level off, due to 
the fixed lowest charge value Z=3 (fig. [17]). Following models, this indicates 
that, as expected, the thermal energy increases with the available energy. For 
better evidencing the evolution with the system mass, several cuts at given 
available energies are performed in fig. [THl one observes that the largest frag- 
ment charge linearly rises with the system mass up to mass 180-190 above 
which a much smoother variation is observed (1-2 units for 100 masses). There- 
fore the independence of the Z^ax distributions on the system only holds for 
heavy enough systems, and starting at higher mass (or charge) than in the 
MMMC calculation. 



Zmax also presents some specific dynamical properties. As shown in [78|, llOl 



for selected compact multifragmenting sources, its average kinetic energy is 
smaller than that of other fragments with the same charge. The effect was ob- 
served whatever the fragment multiplicity for Xe-|-Sn between 32 and 50 MeV/nucleon 
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and for Gd+U at 36 MeV/nucleon. The fragment-fragment correlation func- 
tions are also different when one of the two fragments is Z^^x (see section [7]). 
This observation was connected to the event topology in multifragmentation, 
the heavier fragments being systematically closer to the centre of mass than 
the others. The BOB stochastic mean field simulations, which reproduce a 
large number of the static and dynamical properties of multifragmentation, 
confirm this interpretation, as shown later in fig. [29l It will be shown in the 
next sections that Z^^x also displays scaling properties or bimodal distribu- 
tions, indicating that it may be considered as the order parameter of a phase 
transition in nuclei. 



J^-.5 Summary on size effects 



What has been observed about the influence of the size of the fragmenting 
system on static variables? Firstly the maximum fragment (Z>3) multiplicity 
is proportional to the system mass, for central collisions as well as for QP's. 
In central collisions, this maximum occurs at an energy increasing with Asys 
up to Asys ~ 100-150; the scaling of Mf with Zbound is verified only for very 
heavy systems, and not for a light system with Zsys=56; in that case it can be 
stressed that a very large fraction of Zf,ound is exhausted by He nuclei (more 
than 40% in the explored range). The charge of the largest fragment, for a 
given available energy, increases with Asys up to Agys ~ 190. Thus there is not 
a single system mass which would be a milestone for all the variables. Finally, 
because of its special properties, the largest fragment must be distinguished 
from the others in the partition; being sure of the charge completeness of the 
detected events is mandatory in that aim. 



5 Theoretical descriptions of nuclear fragmentation and compari- 
son to data 



Among the existing models some are relate d to statistica l descriptions based 
on multi - body pha se space calculations [37|, 1 1021 . Il03l . [lOJ, llOSl . Il06l . llOTl . Il08 . 



109, 



111 



112l | whereas others describe the dynamic evolut ion of systems 



resu l ting from colli s ions betw een nuclei via molecular dynamics [113 



124 



114, 



TTel . [m , Jnsl, fiii . [i2ol . Il2ll | or stochastic mean field approaches |122| . 



115, 



123, 



m Im [13 [13 Ea H, [ml, [Tsi, [TsJ- The first approach uses the 



techniques of equilibrium statistical mechanics (counting of microstates) with 
a freeze-out scenario and has to do with a thermodynamical description of 
the phase transition for finite nuclear systems. The second, in principle more 
ambitious, completely describes the time evolution of the collisions and thus 
helps in learning about nuclear matter (stiffness of the effective interaction 
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and in- medium NN cross-sections), its phase diagram, finite size effects and 
the dynamics of the phase transition. Therefore it is highly instructive to 
compare results of the two types of models to experimental data. 



5.1 Statistical ensembles and models 



Three different ensembles are used in statistical models; rigorously only the 
microcanonical ensemble is adapted to describe isolated excited finite nuclei. 
The grandcanonical or macrocanonical ensemble corresponds to the rougher 
description where the system can exchange particles as well as energy with 
a reservoir. In this ensemble the temperature and the chemical potential are 
fixed variables. The total energy but also the total number of nucleons and 
the total charge fluctuate from channel to channel. Constraints are only on 
the average mass and charge of the system. This ensemble is generally used 
for infinite systems and in relativistic quantum systems where particles are 
created and destroyed. However it is a good approximation for hot nuclei when 
one is not interested in an event-by-event analysis and only wants to calculate 
mean values at very high excitation energies (>6-7 Me V/ nucleon) where the 
number of particles associated to deexcitation is large |33, [lo3, H, ffl, [TToj 



(see subsection 12. 2p . A second ensemble, the canonical ensemble, is used to 
describe a system with a fixe d number of particles in contact with a heat 



reservoir at fixed temperature [l06l . Il34j . Il35l | . Here the total energy fluctuates 



from partition to partition and only the mean value of the total energy is 
fixed. Considering hot nuclear systems and using the Fermi gas model, one can 
estimate the standard deviation of the excitation energy distribution, 
function of the temperature T and of the number of nucleons A. It is given by 
the relation a = AE* / \/ AT which shows that the canonical ensemble becomes 
a reasonable approximation for A>200 and temperature T>6 MeV. Within 
the formalisms of both these ensembles many studies have been performed to 
derive mean properties and to discuss phase transitions in terms of intensive 
variables such as temperature. The third ensemble, the microcanonical one, 
is the most relevant for studying isolated systems like nuclei. It i s used to 
describe a system which has fixed total energy and particle number 104l . Il05 . 

In this ensemble the temperature is no longer a natural concept and 



a microcanonical temperature can be introduced through the thermodynamic 
relation: T~l^j,^ = dS/dE. Such an ensemble is fully appropriate if one wants 
to study, for example, partial energy fluctuations (see section [9]) and/or to 
perform analyses on an event-by-event basis. Results have to be discussed 
as a mixing of microcanonical ensembles in order to be compared to those 
of canonical ensembles. Numerical realizations are possible after elaborating 
specific algorithms based on the Monte Carlo method. Finally one can conclude 
about the choice of the different ensembles by saying that the excitation energy 
domain, the size of the system, the pertinent observable to study and the 
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event sorting chosen impose (or permit with some approximation the statistical 
ensemble to be used. For comparison with data additional co nstraints (volume, 
pressure, average volume. . . ) are added to these ensembles [1361 . 137 |. 



5.2 Statistical descriptions of multifragmentation 



A statistical theory of multifragmentation is supposed to predict partition 
probabilities at statistical equilibrium. Thus the weight of a given break- 
up channel /, i.e. the number of microstates leading to this partition, is 
determined by its entropy, Ar/=expS'/, within the microcanonical frame- 
work. Statistical model event generators have been developed for comparisons 
with experiments and in such models the set of fragments corresponding to a 
given partition is considered to be distributed randomly in a freeze-out vol- 
ume equal to 3-10 times the corresponding volume at normal density; the 
freeze-out stage is a ssumed to follow a compression-expansion or/and a ther- 
mal expansion |l38l] phase. These fragments interact via Coulomb forces and 
are endowed with some internal excitation energies and initial thermal ve- 
locities. Radial expansion velocities, fully decoupled from thermal properties, 
are also added in some models. The subsequent evolution including sequen- 
tial decays of primary fragments is performed, preserving in some models 
space-time correlations. Finally generated events are filtered to account for 
the experimen tal device . Detailed presentations of models can be found in 



references [17|, llOSl . IllOl . Illll . Ill2| . Within such an approach the input pa- 
rameters: mass and charge of the system at break-up density, its excitation 
energy, its volume at freeze-out and the eventual added radial expansion have 
to be backtraced to experimental data, estimated from dynamical simulations 
or derived from data related to properties of systems at break-up. Statistical 
models are very useful to compare data sets and to produce comparisons or 
systematics. The following examples illustrate their different utilizations. 
In fig. [19] the fragment charge distribution me asure d for central ^'^^Xe+^°'^Cu 
collisions at 30 MeV/nucleon incident energy 139l | is compared to the Berlin 
Multifragmentation Model (also called MMMC-Microcanonic al Me tropolis Monte 
Carlo) [loi] and to the sequential statistical model GEMINI [m^. The BMM 
calculations give good qualitative agreement with the experimental charge dis- 
tribution over the entire Z range measured. Input parameters of the source 
(A=177, Z=76 and 750 MeV excitation energy) were estimated from dynam- 
ical simulations at a freeze-out time corresponding to a volume of 5.6 Vq 
(i?i?o=ll-8 fm for the sphere including all the products at freeze-out). Re- 
sults with a larger freeze-out volume (10.7 Vq) which show a better overall 
agreement are also shown. Conversely, the GEMINI calculation underpredicts 
the yields of low Z by an order of magnitude and predicts a nearly flat charge 
distribution. A comparison with the most popular model, namely the Statis- 



tical Multifragmentation Model SMM 110| . is shown in fig. [20] for the largest 
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fragment and the fragment multiplicity as a function of the thermal excita- 
tion energy of the source. It concerns multifragmentation of 1 GeV/nucleon 
Kr, La and Au projectiles pro duced in collisions with a carbon target, stud- 
ied by the EOS collaboration 14ll . Il42( |. Inputs for SMM (thermal excitation 
energy, mass, charge and freeze-out volume of the source) were derived from 
data including p-p correlations for freeze-out volume and the ISABEL cas- 
cade model for estimate of preequilibrium neutrons, which were not detected. 
The ranges of mass and charge of the SMM sources which cover the explored 
thermal excitation energy domains are A=80-30 and Z=34-16 for Kr, A=130- 
80 and Z=54-26 for La and A= 190- 120 and Z=76-48 for Au. The standard 
value of SMM (16 MeV) was used for the inverse level density parameter 
and derived freeze-out volumes correspond to three times the volume of Kr, 
La and Au nuclei at normal density. In SMM the definition of the freeze-out 
volume differs from that of BMM and includes only the centre of all the frag- 
ments. Very good agreements are observed between data and SMM for the 
largest fragment over the whole range of thermal excitation energy whereas 
fragment multiplicities are overestimated by about 20% in the excitation en- 
ergy range 5-10 MeV/nucleon. Finally fig. [2T] shows different static and dy- 
namic fragment observable distributions measured for central Xe-|-Sn collisions 
at 32 MeV/nucleon incident energy and compared again with SMM. To get 
this agreement the input parameters of the source are the following: A=202, 
Z=85 as compared to A=248 and Z=104 for the total system, freeze-out vol- 
ume 3Vo, partitions fixed at thermal excitation energy = 5 MeV/nucleon and 
added radial expansion energy of 0.6 MeV/nucleon. As we have shown with 
those examples, but it is a general trend, the observed properties of fragments 
are compatible with the hypothesis of sources in thermal equilibrium which 
undergo multifragmentation. 



5.3 Dynamical descriptions of multifragmentation 



Beside statistical descriptions, there are microscopic frameworks that directly 
treat the dynamics of colliding systems such as the family of semi-classical 
simulations based on the nuclear Boltzmann equation (the Vlasov-Uehling- 
Uhlenbeck (VUU), Landau- Vlasov (LV), Boltz mann - Uehl i ng-Uhlen beck (BUU) 
or Boltzmann-Nordheim- Vla sov ( BNV) co des 



molecular d ynamics ( CMP ) |12ll . Il48l . Il49l . Il50l |. quantum molecul ar d ynam- 



143 . [lisl . \im [mJ), classical 



ics (QMD) jll3 . 1141 . 115 |. fermionic mol ecular dynami cs (FMD) jll7 |. anti- 
symmetrized molecular dynamics (AMD) jll6 . 118, 119 1 and stochastic mean 
field approaches related to simulations of the Boltzmann-Langevin equation 



122l . Il23l . Il24l . Il26l . Il27l . Il28| . Boltzmann type simulations follow the time evo- 



lution of the one body density. Neglecting higher than binary correlations, they 
ignore fluctuations about the main trajectory of the system (deterministic de- 
scription), which becomes a severe drawback if one wants to describe processes 
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Fig. 19. Experimental charge distribution (full points) in central events compared 
with unfiltered (dot-dashed curves) and filtered (solid curves) predictions by the 
BMM model wi th t wo freeze-out radii (top and central panels, see text) and by 
GEMINI (from fl39(). 
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Fig. 20. Left panel: size of the largest fragment as a function of thermal excitation 
energy, E^j^ for multifragmentation of 1 GeV/nucleon Kr, La and Au (data) and 
SMM. Right panel: idem for average fragment multiplicities (from 142| ). 



involving instabilities, bifurcations or chaos expected to occur during the mul- 
tifragmentation process. Such approaches are only appropriate during the first 
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Fig. 21. Comparison of experimen- 
tal data (Xe+Sn compact single 
sources produced in central colli- 
sions-32 MeV/nucleon) with SMM 
simulations. The lines are for data 
and symbols for SMM (all fragments, 
except open circles and dashed lines 
which refer to the largest fragment of 
each partition). Zbound represents the 
sum of the charges of all fragments 
(from 0)- 

stages of nuclear collisions, when the system is hot and possibly compressed 
and then expands to reach a uniform low density. They become inadequate 
to correctly treat the fragment formation, and for the description of multi- 
fragmentation it is essential to include fluctuations. This is done in quantum 
molecular dynamics methods and in stochastic mean field approaches. 



Fig. 22. Fragment (left panel) and 
heaviest fragment (right panel) 
charge distributions for central 
^^'^ Au+^^'^ Au collisions at different 
incident energies. Points refer to 
data and full lines to fragment from 
QMD events. Fragments from QMD 
events filtered by the experimental 
device cor respond to dashed lines 
(from |l4i|). 



5.3.1 Quantum molecular dynamics: QMD and AMD simulations 

QMD is essentially a quantal extension of the molecular dynamics approach 
widely used in chemistry and astrophysics. Starting from the n-body Schrodinger 
equation, the time evolution equation for the Wigner transform of the n-body 
density matrix is derived. Several approximations are made. QMD employs 
a product state of single-particle states where only the mean positions and 
momenta are time-dependent. The width is fixed and is the same for all wave 
packets. The resulting equations of motion are classical. Also the interpretation 
of mean position and momenta is purely classical and the particles are consid- 
ered distinguishable; this simplifies the collision term which acts as a random 
force. All QMD versions use a collision term with Pauli blocking in addition to 
the classical dynamics. Some versions consider spin and isospin and others do 
not distinguish between protons and neutrons (all nucleons carry an average 
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charge). Finally one has to stress that for most of the QMD versions a statisti- 
cal decay code mus t be coupled to describe the long time dynamics. However 
for the code of ref. [1141] there is no ne ed t o supplement the QMD calculations 



by an additional evaporation model |15l| . This code was very recently used 



in a rather complete comparison with data measured for central ^ ^'^Au +'^^'^Au 



collisions over the incident energy range 60-150 MeV/nucleon |143| . Using 



a soft equation of state (incompressibility Koo= 200 MeV), static properties 
like fragment multiplicity and charge distributions are rather well reproduced, 
particularly for the higher incident energies where QMD codes are certainly 
better adapted. Figure [22] illustrates the comparison for fragment and heaviest 
fragment charge distributions. 

An antisymmetrized version of molecular dynamics (AMD) was constructed 
by incorporating two-nucleon collision process as the residual interaction into 
the fermionic molecular dynamics (FMD). AMD describes the system with 
a Slater determinant of Gaussian wave packets and therefore can describe 
quantum-mechanical features. However, in the dynamics of nuclear reactions, 
there may be other phenomena caused by the wave packet tail that are com- 
pletely lost in AMD due to the restriction of the single-particle states. So an 
improvement was realized (called AMD-V) with the stochastic incorporation 
of the diffusion and the deformation of wave packets which is calculated b y th e 



Vlasov equation without any restriction on the one-body distribution [118 
More recently the quantum branching process due to the wave packet diffusion 
effect was treated as a random term in a Langevin-type equation of motion 
whose numerical treatment is much easier. Moreover a new approximation 
formula was also introduced in order to evaluate the Hamiltonian in the equa- 
tion of motion with much less computation t ime than the exact calculation, so 



that systems like Au-|-Au are now treatable |152l |. As for QMD the stiffness of 
the effective interaction and the in-medium NN cross-section are both impor- 
tant ingredients for determining the degree of agreement with experimental 
data. In order to test the sensitivity of the ingredients, a detailed study of 



reaction dynamics and multifragmentation was done in ref [81[ by comparing 
AMD-V calculations with data from heavy-ion reactions around the Fermi 
energy. Figure [23] presents multiplicity distributions of selected particles and 
fragments produced in central ®'^Zn+^^Mo collisions at 47 MeV/nucleon in- 
cident energy. Thin solid, dashed and thick solid lines indicate the results of 
Soft-|-A^A^emp, SiiS+NNemp and SUS+NNlm respectively. As previously ob- 
served with QMD a better global description of data is obtained with a soft 
equation of state (incompressibility Koo= 228 MeV). We refer the reader to 
ref. ISli for more information on in-medium NN cross-section. 



5.3.2 Stochastic Mean Field approach: Brownian One-body Dynamics (BOB) 

In many domains of physics a diffusive behaviour is described by transport 
theories which were originally developed for Brownian motion. The effects of 
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Fig. 23. Multiplicity distributions of 
selected particles and fragments pro- 
duced in central ^^Zn+^^Mo collisions at 
47 MeV/nucleon incident energy. Exper- 
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Fig. 24. Comparison of experi- 
mental data (central collisions: 
Gd+U-36 MeV/nucleon and 
Xe+Sn-32 MeV/nucleon) with BOB 
simulations for charge and multiplicity 
distributions of fragments (top and 
middle panels) and for their average 
kinetic energies (bottom panel). The 
symbols are for data and the lines for 
BOB simulations. Light grey lines and 
triangles stand for Gd+U and black 
lines and circles for Xe+Sn (adapted 
from [91]). 



the disregarded degrees of freedom are simulated by a random term in the 
dynamics of the retained variables. That idea is the starting point of the so- 
called Boltzmann-Langevin equation (BLE): df/dt = ih\f], f} + I\.f] + SI \f] 
which was introduced for heavy-ion collisions in references jl22 . 123 . 124 1. / 
is the one-body phase space density. The first term on the r.h.s. produces the 
collisionless propagation of / in the self-consistent one-body field described 
by the effective Hamiltonian. The second term, called collision term, repre- 
sents the average effect of the residual Pauli-suppressed two-body collisions; 
this is the term included in LV, BUU and BNV simulations. The third term 
is the Langevin term which accounts for the fluctuating part of the two- 
body collisions. Exact numerical solutions of the BLE are very difficult to 
obtain and have only been calculated for schematic cases in one or two dimen- 
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sions [125j. Therefore various approx imat e treatments of the BLE have been 
developed. The basic idea of BOB [l29'| is to replace the fluctuating term 
by SI[f] = —6F[f].df /dp where 5F(r, t) is the associated Brownian force 
(< 6F >= 0). Since the resulting Brownian one-body dynamics mimics the 
BL evolution, the stochastic force is assumed to be local in space and time. 
The strength of the force is adjusted to reproduce the growth of the most un- 
stable modes for infinite nuclear matter in the spinodal region (see section [8]). 
Quantal fluctuations connect ed w ith collisional memory effects are also taken 



into account as calculated in 153 



An extensive comparison data-BOB was made for two very heavy fused sys- 
tems produced in Xe-|-Sn and Gd+U central collisions which undergo multi- 
fragmentation with about the same excitation energy (~ 7 MeV/nucleon) [t^, 
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91j. Stochastic mean-field simulations were performed for head-on colli- 
sions with a self-consistent mean field potential chosen to give a soft equation 
of state (Koo= 200 MeV). The finite range of the nuclear interaction was taken 
into account using a convolution with a Gaussian function with a width of 0.9 
fm. A term proportional to Ap in the mean-field potential was added; it allows 
to well reproduce the surface energy of ground-state nuclei, which is essential 
in order to correctly describe the expansion dynamics of the fused system. 
In the collision term a constant NN cross-section value of 41 mb, without in- 
medium, energy, isospin or angle dependence was used. As a second step the 
spatial configuration of the primary fragments, with all their characteristics 
as given by BOB, was taken as input in a statistical code to follow the frag- 
ment deexcitation while preserving space-time correlations. Finally the events 
were filtered to account for the experimental device. These simulations well 
reproduce the observed charge and multiplicity distributions of fragments (see 
fig. [21]). Particularly the independence of the charge distribution against the 
mass of the system experimentally observed was recovered 76|. More detailed 
comparisons of the char ge d istributions of the three heaviest fragments also 



show a good agreement [91[. Kinetic properties of fragments are rather well 
reproduced for the Gd+U system, whereas for Xe+Sn the calculated energies 
fall ~ 20% below the measured values. We also refer the reader to sections [7] 
and [8] for further comparisons. 



5.3.3 Concluding remarks on dynamical simulations 

To conclude on dynamical descriptions of multifragmentation one can make 
a few general comments. Better agreements between data and calculations 
are generally observed with a soft equation of state (incompressibility Kqo ~ 
200-230 MeV). One can however not trust the in- medium NN cross-sections 
extracted as the self-consistency between the mean field potential and the 
two-body collision term is not fulfilled in models. QMD and AMD calcula- 
tions lead to too much transparency at low incident energies (typically below 
50-100 MeV/nucleon). For AMD-V and BOB calculations at incident ener- 
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gies around 35 MeV/nucleon a maximum density of 1.2-1.4 po is observed at 
30-40 fm/c after the beginning of central heavy- ion colhsions and the normal 
density is recovered around 70 fm/c. Thermal equilibrium times are found in 
the range 100-140 fm/c after the beginning of collisions, well before freeze-out 
configurations (200-300 fm/c). Primary fragments exhibit an equal or almost 
equal excitation energy per nucleon of 3-4 MeV. The mechanism of fragment 
production differs depending on model type. In molecular dynamics models 
fragments are preformed at early stages close to the normal nuclear density 
whereas in stochastic mean field calculations, as BOB, fragment formation 
is linked to the spinodal decomposition mechanism: mononuclear systems at 
low density (~0.4po) formed at around 100 fm/c develop density fluctuations 
during about 100 fm/c to form fragments. More constrained observables re- 
lated to the formation of fragments by spinodal instabilities will be discussed 
in section [SI 

Last point, dynamical calculations exhibit radial collective energies for frag- 
ments with average values in the range 0.1-2.0 MeV/nucleon for heavy-ion 
collisions in the Fermi energy domain which fairly agree with values derived 
from experiments (see subsection 14.31) . 



5.4 The link between dynamical and statistical descriptions 



While both descriptions show reasonable agreement with data in reproducing 
average static and kinematic properties of fragments, sharp conclusions on 
multifragmentation scenarii can not yet be derived. However from dynamical 
descriptions like AMD-V and BOB we learnt that, in average, thermal equilib- 
rium was reached before the breaking stage, which is also the hypothesis of sta- 
tistical approaches. A tentative to derive more precisely at what tim e, durin g 



the reaction, the statistical description takes place was done in refs. [78|, 1154 
in BOB simulations of Xe+Sn central collisions at 32 MeV/nucleon incident 
energy, the volume of the system was calculated every 20 fm/c from 100 to 
250 fm/c to be compared with the freeze-out volume of 3Vo as input in the 



SMM code |78l. Il54l |. The result is the following: at 200 fm/c, during the for- 
mation of fragments, a freeze-out volume of ^3Vq (calculated a la SMM) is 
reached and the characteristics in size of the source are A=190 and Z=80, 
in good agreement, within 10%, with SMM inputs. Very recently a remark- 
able result was also obtained by checking the consistency of predictions be- 
tween the dynamical evolution of Xe+Sn central collisions at 32 MeV/nucleon 
via a stochastic mea n fie ld approach and a microcanonical multifragmenta- 
tion model (MMM) jl55 |: a statistically equilibrated stage was identified at 
~140 fm/c just intermediate between the beginning of the spinodal decompo- 
sition and the freeze-out configuration. In both cases a rather coherent link 
between dynamical descriptions, like stochastic mean field approaches, and 
statistical descriptions is derived. 
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6 Calorimetry and thermometry 



The knowledge of the excitation energy and temperature, as weU as of its 
numbers of neutrons and protons (mass and charge) is mandatory for any 
description of a system in terms of thermodynamical variables. Determining 
these quantities remains the most challenging task for the groups involved in 
the study of multifragmentation. 

As previously mentioned, in nuclear collisions the formed multifragmenting 
systems are accompanied by a quite abundant preequilibrium emission. These 
early emitted products should not be included in the calculation of the mass 
and of the excitation energy of the fragmenting system; it is unfortunately 
difficult to unambiguously attribute the observed final products to preequi- 
librium or to the fragmenting system. This is particularly true for neutrons 
- when measured - and for light charged products, H and He and to a lesser 
extent Li and Be isotopes. Conversely the properties of the heavier fragments 
indicate that they do originate from the fragmenting system. Preequilibrium 
is removed either with the help of models, or through angular and energetic 



properties of the observed products [77|, l91 



6. 1 Calorimetry 



All procedures for obtaining the excitation energy of a fragmenting source, ob- 
served with a 47r array, are based on the determination of its velocity. For cen- 
tral collisions the reaction centre of mass velocity is most often chosen whereas 
the quasi-projectile velocities are identified to either that of the biggest frag- 
ment, or that of the subsystem containing all the fragments (Z>3 or 5), for- 
ward emitted in the centre of mass. The excitation energy, E*, of the source is 
then calculated event by event with the relation E* = J^M^p Ecp + J2m„ En — Q. 
Ecp and En are respectively the kinetic energies of charged products and neu- 
trons belonging to the source, Q is the mass difference between the source 
and all final products. Energies are expressed in the source reference frame. 
Mcp is in most cases the detected multiplicity of charged products. The energy 
removed by gamma rays is small and most often neglected in the calculation. 
The decision of including or not one observed charged product in the source 
differs with the experimental apparatus and the type of collision under study. 

CC for central heavy-ion collisions, all fragments with Z>3 (or 5) are attributed 
to the source. Preequilibrium in that case is mostly forward/backward emit- 
ted, and indeed the angular distributions of the light products appear isotropic 
between 60 and 120°. The charge, mass and energy contributions of these 
particles are doubled for the calculation of the characteristics of the source 
(i). Another possibility, to account for the detector inefficiency, is to calcu- 
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late the charge, mass and energy of the anisotropic part, and to remove it 
from those of the composite system (ii). 
QP for quasi-projectiles the most important contamination comes from mid- 
rapidity products and several techniques are used for the QP reconstruction, 
i) All fragments forward emitted in the reaction cm. system are attributed 
to the QP. Vari a nts c onsist either in putting a low velocity cut for the lighter 



fragments |156l . Il57l |. or in kee ping only events with a compact fragment 
configuration in velocity space 74, ISq . Then twice the light elements in 
the QP forward hemisphere are added, ii) Fragments are treated as above, 
but particles are attributed a probability to come from the QP emission. 



either using a 3-so urce fit [42j, or by taking a well characterized subspace as 



reference |50l . Il59| . The velocity of the QP is then recalculated by including 



all its components. 

hIC finally in hadron-induced collisions, products emitted from the source are 
chosen from energetic considerations, by excl uding those with an ener gy p er 



nucleon above a given threshold either fixed |160| or varying with Z |161 



All those procedures assume forward-backward symmetry of particle emission 
in the source frame. For QP's the symmetry of the source emission may be 
questionable wh en h i ghly excited QP's and QT's start emitting right after 
their separation 1621 . Il63| : the close proximity of the partner deforms phase 
space and emission is favoured between the QP and QT. This possible effect 
is generally ignored. 

Once counted the charged products, the charge of the source is known. A first 
uncertainty is introduced in calculating the associated mass, as that of heavy 
fragments is not measured. A single mass is attributed to all nuclei with a 
given atomic number, either that of the most stable species, or th at derived 
from formulae existing in the literature (EPAX |l64 | or EAL 165| ). At that 
point neutrons must be included. Except in experiments using a neutron ball, 
neither their multiplicity not their energy is known. The neutron number is 
thus calculated by assuming that the source has the same N/Z ratio as the 
total system (central or hadron-induced collisions) or as the projectile. The 
average neutron energy is taken equal to the average proton energy over the 
event sample, removing some Coulomb barrier. Note that with neutron balls 
only the neutron multiplicity is measured, at the price of a poor geometrical 
coverage for charged products. In that case corrections accounting for the un- 
detected particles and neutrons are made 
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In central heavy-ion reactions, the excitation energy domain populated is 
narrow : cte* ~0.7-1.25 MeV/nucleon; the width includes experimental effects 
(detector efficiency and resolution), calculation assumptions and physical ef- 
fects (pre-equilibrium) . Conversely, in hadron-induced reactions as well as in 
quasi-projectile studies, a broad domain of excitation energy is populated, 
proportionally to the partial cross section, function of the impact parameter. 
However, due to on-line trigger effects, very low energies are poorly sampled, 
by particular events matching the trigger requirement; indeed neutron emis- 
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Fig. 25. Top : Convoluted, unconvoluted and 
experimental excitation energy distributions ob- 
tained in vr+Au reactions. Bottom : Ratio of un- 
convoluted-to-convoluted distribution func- 
tion o f the excitation energy per nucleon. Adapted 
from 166l | 



Fig. 26. Relative share of 
excitation energy for vari- 
ous components of the re- 
construction procedure as a 
function of the excitatio n en - 
ergy per nucleon. From [16d | 



sion - not detected - is dominant in this region. At the other end of the dis- 
tribution, the very high energies probably result from significant fluctuations. 
In all cases the reliable domain extends from about 2 to 8 MeV/nucleon: for 
example in fig. [25] obtained in vr+Au reactions, the excitation energy distri- 
bution is unconvoluted assuming Gaussian fluctuations. More than half of the 
1% of events abov e th e vertical dotted line have an energy overestimated by 



1-2 MeV/nucleon [166 



How reliable are the energies so obtained ? Because of compensation of the 
errors on the mass and on the energy, the energy per nucleon is a more robust 
experimental observable than the energy. By comparing values obtained by 
different metho ds fo r quasi-projectiles, differences on E* /A smaller than 10% 



were found |74l. Il67l |. From simulations with the HIPSE event generator, the 
reconstructed values were found to differ from the true values by less than 
10%, e xcep t for very peripheral collisions where the discrepancies are much 



larger [159[. In central collisions, excitation energies slightly smaller than the 
available energies are generally found, which is what can reasonably be ex- 
pected. It was verified in the INDRA Xe+Sn data that procedures CCi and 
CCii give the same results when a high degree of completeness is required 
for CCi, e.g. 90% of the system charge. For lesser completeness (80%), the 
difference between both types of calculation increases with the incident en- 
ergy, reaching 1 MeV/nucleon (10%) at 50 MeV/nucleon. The main source 
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of uncertainty in the calculation of E* comes from the neutron terms. How- 
ever, compensation occurs in the calorimetry equation between the kinetic and 
the mass balance terms; indeed the weight of these two terms is simila r for 



quasi-projectiles and in hadron-induced reactions {Q/E* ~30-36% [7J, [166 
- see fig. |26|) : in central collisions the Q term accounts only for ~20% of the 
excitation energy. 



6.2 Temperature measurements 



Two recent rev iews exte nsively describe the methods used for temperature 



measurements 



168l . Il69l . and references therein]. A brief summary will be 
given here. The concept of temperature for a nucleus, which is a microscopic, 
isolated Fermionic charged system, is not a priori obvious. According to sta- 
tistical mechanics, temperature can be defined as = dS{E, N,V)/dE, 
where S, E, V, N, are the entropy, energy, volume and (fixed) number of 
particles of the system. This definition is valid only if the system is in statis- 
tical equilibrium and its density of states is known. For compound nuclei at 
low excitation energy, both conditions are fulfilled, which might not be the 
case at higher energies, and for fragments resulting from multifragmentation. 
Moreover no probe can be used to measure the temperature of these small 
systems, it has to be derived from the properties of particles that they emit 
during their cooling phase. The abovementioned issues of identifying an equi- 
librated emitting source by distinguishing pre-equilibrium particles from those 
emitted by that source is therefore common to all the methods. Three families 
of methods are used to "measure" temperatures. 

(1) Kinetic approaches. Historically, temperatures of compound nuclei were 
derived from the slopes of the kinetic energy spectra of the emitted neu- 
trons or charged particles that they evaporate, as the spectra can be 
fitted with Maxwell-Boltzmann distributions. At higher energies, when 
long chains of particles are emitted, the obtained result is an average 
over the deexcitation chain, and may differ from one particle to another, 
depending on the emission sequence. To retrieve the initial temperature, 
it was proposed to subtract from the spectra those of p articl es coming 
from the same nucleus formed at lower excitation energies For mul- 



tifragmenting systems, the slopes of light product spectra lead to very 
high "temperatures", and do not probably reflect only the thermal prop- 
erties of the system, but also the collective energies coming from the 
dynamics of the nuclear collision. It was recently proposed to derive the 
temperat ure of the fragmenting system from the slope of the thermal hard 
photons flT?], which have the advantage of being insensitive to Coulomb 
field and final-state effects: temperatures close to 7 MeV are for example 
obtained for central Xe-|-Sn collisions at 50 MeV/nucleon. 
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(2) Population of excited states. The underlying idea for this method is that 
the population of the excited states of a system in statistical equilibrium 
is given by the temperature of the system and the energy spacing, AE, 
between the levels. This definition in itself bears the limits of the method: 
when the temperature is higher than AE, the ratio between the popula- 
tion of two states saturates. Considering particle-unbound states is thus 
interesting as it allows to measure higher temperatures, and the popula- 
tion ratio should in that case be less influenced by secondary decays. 

(3) Double ratios of isotopic yields. This method uses the yields of different 
light isotopes produced by the system. It was developed in the grand- 
canonical approach, and is valid for systems at densities low enough 
to make fragment nuclea r int eraction negligible, thus the composition 



of the system is frozen |172l ]. The basic assumption is that free nu- 
cleons and fragments are in thermal equilibrium within an interaction 
volume V. The density of an isotope reads: p{A, Z) = N{A, Z)/V = 
A^/^u;{A, Z)XtI exp(/i(^, Z)/T), where uj is the internal partition func- 
tion of particle(y4, Z), /i its chemical potential and A the thermal nucleon 
wave length. The condition of chemical equilibrium allows to define the 
chemical potential of a species in terms of those of free neutrons and 
protons and of its binding energy. Using two sets of two nuclei differing 
only by one nucleon, the temperature is derived from the double yield 
ratio, the binding energy differences and the partition functions only, the 
other terms disappear. The results depend on the validity of the assump- 
tions: is the grandcanonical ensemble relevant in the studied case? Is the 
system really in thermal and chemical equilibrium? The considered par- 
ticles should be present at freeze-out, and not produced by secondary 
decays. Another problem lies in the calculation of the binding energies 
which might depend on density and temperatur e. D ifferent corrections 
were proposed to account for finite-size effects |l73| or secondary de- 



cays 



(4) Isospin thermometer approach. At relativistic energies, it was proposed 
to derive the temperature at freeze-out from the isotopic distributions 
of the final residues 176| . The assumptions are that the fragments at 
freeze-out have the same N/Z ratio as the projectile; their thermal en- 
ergy is dissipated through an evaporation cascade, which is reconstructed 
with the help of evaporation codes. The results thus entirely rely on the 
reliability of these codes. 



A recent study on the reliability of temperature measurements for Xe QPs 
from Xe+Sn collisions at various energies shows that, whatever the method 



used, temperatures are determined at best within 10-20% |l59 
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Fig. 27. Caloric curves for three selected regions of mass. The initial temper atur es 
are plotted vs the thermal excitation energies of the systems. (Adapted from [l77|]). 

6.3 Caloric curves 



The combined and independent determinations of the thermal excitation en- 
ergy and of the temperature of a system allow to draw the caloric curve, which 
gives the rela tion between these variables. The first example of caloric curve 



was given in |178l | and taken as the evidence of the occurrence of a liqu id-gas 



type phase transition. Many data followed, which were compiled in |177l | , from 
which is extracted fig. [23 Care was taken to recover the initial temperature 
of the systems by applying corrections, as mentioned in sect. 16.21 to account 
for the deexcitation cascade, or with the help of a statistical model. The pub- 
lished excitation energies were considered as pure thermal energy, except in 
the case of a heavy ion reaction for which the thermal part was evaluated 
with the SMM model. The general trend of these curves, classified in different 
mass zones, is first the well-known rising part attributed to the Fermi gas 
{E* = aT"^, with a the level density parameter), followed by a plateau. The 
temperature at the plateau, as well as the energy at which it is first reached, 
were seen to decrease for larger masses. This reminded of the evolution of the 
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predicted limiting temperatures caused by Coulomb instabilities |179l . II8O 
More recently, caloric curves were publis hed which use the temperatures ob- 
tained from hard thermal photon spectra [l71 . 181 |. The results do not clearly 
exhibit a plateau, but above 3-4 MeV/nucleon of excitation the caloric curves 
fall below that expected for the Fermi gas. Finally a theoretical study in the 
microcanonical framework of the MMM model calib rated the different isotopic 
temperatures against the true microcanonical one 173| : the authors found a 
universal relation for all masses, and used it to re-evaluate some published 
caloric curves. The corrected curves all exhibit three parts, the Fermi gas one, 
a more or less broad plateau followed by a linearly rising part attributed to a 
"classical" gas region. Further discussion on caloric curves will be presented 
in sect. [8] in relation with phase transitions. 



7 Freeze-out properties 



The concept of freeze-out was first introduced as a starting point in statis- 
tical models. It can be defined as a configuration for which mutual nuclear 
interactions between primary products of a multifragmenting source become 
negligible. Moreover, for the sake of simplicity, those products are supposed 
to be spherical and having eventually recovered normal density. It was shown 
recently, using dynamical simulations as BOB, that such a simple geometri- 
cal picture can be relevant on the event by event basis |182i] . The freeze-out 
configuration characterizes the new physical state produced by the excited fi- 
nite system. In particular the freeze-out volume appears as a key quantity and 
its knowledge is particularly important in the extraction of thermostatistics 
observables as the microcanonical heat capacity (see section [H]). Relative ve- 
locity correlations including fragments and particles as well as different types 
of simulations compared to data have been used to progress with the difficult 
task of determining the freeze-out volume. 



7.1 Fragment velocity correlations and event topology at freeze- out 



Fragment-fragment relative velocity correlations functions were widely used 
in the 90 's to derive information on emission time scales and the disassem- 
bling source volume. Most of t he pa pers published on the subject make use 



of the formalism developed in [183| . where arguments are presented which 
justify a classical treatment of such correlations. The importance of Coulomb 
effects and momentum conservation laws is underlined. In this picture, it is 
suggested to mix correlations obtained with fragments of different sizes by 
replacing the relative velocity, Vrei, between two fragments of charge Zj and 
Zj by a reduced velocity, Vred = fre«/ y^-^i + ^jj this formula the mass of 
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the fragments is supposed to be twice their charge, meaning that it is more 
appropriate for hght fragments. The uncorrelated yield necessary to build the 
correlation functio n was gen erally obtained from event mixing between events 



of the same class [18J, |185| . As for particle-particle correlations, one obtains 
information on the space-time extent of the emitting source. A majority of the 
results deal with trajectory calculations, assuming an exponential probability 
distribution of the time delay between emitted fragments, P{t) oc exp—t/r. 
Excitation functions show that, for central collisions, the space-time extent 
decreases when the bombarding energy is raised: times around 300-500 fm/c 
are reported around 30 MeV/nuc leon, while they are between 50 and 100 fm/c 
above 60 MeV/nucleon |l86 . 187]. For multifragmenting quasi-projectiles pro- 
duced in semi-p eriph eral collisions, fragment emission times of 200-500 fm/c 



are reported in |l88 |. for the Cu-|-Au reaction at 50 MeV/nucleon, indepen- 



dently of the assumed volume of the QP source; for Ni QP from the Ni-|-Au 
reaction at 34.5 MeV/nucleon, r decreases from 550 to 75 fm/c when the 
excitati on en ergy increases from 2 to 7 MeV/nucleon and remains constant 
beyond jl89| . Thus on average, QP fragment emission times are comparable 
to those of central collisions, as appears in fig. [3l where the emission times 
are plotted versus the source excitation energy. Finally velocity correlations 
were also studied in hadron-nucleus reactions as a function of the excitation 



energy deposited in the nucleus [98|. Emission times were shown to decrease 
from r=500 fm/c at E*/A=2 MeV to r ~20-50 fm/c for E*/A>5 MeV. The 
time range quoted at high excitation accounts for corresponding associated 
volumes 6-4Vo. Figure [3] shows that the times obtained for these reactions are 
systematically shorter than those measured in nucleus-nucleus collisions, by 



a factor of about 5. The authors of [98| tentatively explain this observation 



by a better source selection in hadron-induced reactions. This would act on 
the excitation energy axis only (through the energy range of included Icp, the 
collective energies). Altogether, these results testify for an increased thermal 
energy deposition in collisions when raising the excitation energy. 
In a few cases, fragment-fragment velocity correlation functions were com- 
pared to a full calculation coupling a dynamical simulation of the collision 
to an after-burner. Phase-space coordinates of the fragments produced in the 
first phase are injected in the second one. Two examples will be given her e, on e 
coupling QMD and SMM for central 100 MeV/nucleon Fe+Au collisions 190|, 
and the other using the stochastic mean field BOB simulation followed by the 
SIMON deexcitation for central 32 MeV/nucleon Xe-|-Sn and 36 MeV/nucleon 
Gd+U reactions 78|. Provided that the models accounts sufficiently well for 
more simple fragment properties, which was verified in both cases, a good 
reproduction of the correlation functions (filtered by the experimental appa- 
ratus) allows to trust the freeze-out volume and the multifragmentation time 
given by the dynamics. For the Fe-t-Au reaction, thermal multifragmentation 
of the quasi-target was shown to be the dominant mechanism. The freeze-out 
volume is given by the SMM parameter k. The overall agreement between 
data and calculations is not very good (fig. [2S]left), nevertheless the rising 
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Fig. 28. Experimental reduced velocity corre lation functions (symbols) for central 
collisions in the r eactions Fe+Au (left US) and Xe+Sn and Gd+U(right [78]) 
for different charges of the fragment pairs. The lines represent calculations from 
QMD+SMM and different freeze-out volumes in the left panel, and from BOB in 
the right panel. See text for details 

part after the Coulomb hole allows to derive freeze-out volumes which in- 
crease between 2 and IOVq (or alternatively shorter emission times) with the 
charge of the fragment pair. Times shorter than 500 fm/c are indicated. That 
small fragme nts are emitted from a hotter and more dense source was also 
suggested in 19l[ | from large angle correlation results. In the case of reactions 
between heavier ions near the Fermi energy, where the formation of compact 
single sources was evidenced, the BOB-I-SIMON model well reproduces the 
three types of pair selections, with no cuts on the fragment energy/velocity: 

i) all fragments considered (Zj j > 5); 

ii) intermediate mass fragments (IMF) 5 < j < 20; 

iii) reduced velocity correlation of the heaviest fragment Z^ax with each of 
the others in the event Z[ < Z^ax- 

Note that Li and Be are not considered in that case. The almost perfect 
matching of the calculated correlation functions with the data, particularly 
for Gd+U (fig. [28] right), led the authors of ref [78| to extract average freeze- 
out volumes (with a definition close to the SMM one) around 4Vo for Xe-|-Sn 
and 8Vo for Gd+U. The freeze-out instant was d efined as that when the 
average fragment multiplicity becomes frozen [isj, namely 200 (240) fm/c 
after the incident partners come to contact for Xe+Sn (Gd+U). The freeze- 
out topologies found in the simulations imply that the heavier fragments are 
located closer to the center, as shown in fig. [211 Note that in ^He induced re- 
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Fig. 29. Fragment positions as a function of the atomic number at the freeze-out 
instant for 32 MeV/nucleon ^■^^Xe+^^^Sn - upper panel, and 36 MeV/nucleon 
^^^Gd+^'^'^U - lower panel. The contour scale is logarithmic. Adapted from (tsI 

actions it was conversely found that the heavier fragment should be randomly 
positioned in a dilute source in order to well ac coun t for the experimental 



correlation functions with a trajectory calculation [191[. This observation may 
sign a difference in the fragmentation process whether it is thermally driven 
as it seems for QP and hadron induced reactions or it follows a compression 
phase as in very heavy-ion central collisions. In that case self similar expansion 
should favour coalescence of primary fragments close to the centre of mass. 



1.2 Fragment-particle correlations and fragment excitation at freeze- out 



Excitation energy measurements of primary fragments which are present at 
the freeze-out stage are also of large interest. They can provide information on 
the degree of equilibration at freeze-out and put strong constraints on freeze- 
out characteristics through a microcanonical description of multifragmentation 
events. Such measurements are derived from fragment-particle correlations. By 
using the correlation technique for the relative velocity between light charged 
particles and fragments, one can extract the multiplicities and average ki- 
netic energies of particles emitted by fragments with a given final charge, 
and then reconstruct the sizes and exci tatio n energies of the primary frag- 
ments. This technique, first proposed in 192| for Xn+Sn central collisio ris at 
50 MeV/nucleon incident energy, was applied from 32 to 5 MeV/nucleon jl93 
and for Kr+Nb central collisions at 45 MeV/nucleon |l94 |. For each bombard- 
ing energy, a constant value of the mean excitation energy per nucleon has 
been found over a primary fragment charge range 5-20, which strongly suggests 
that fragments are on average in thermal equilibrium at freeze-out. This aver- 
age excitation energy, equal to 2.2 MeV/nucleon at the lower incident energy, 
saturates around 2.5-3.5 MeV/nucleon for beam energies 39 MeV/nucleon and 
above. Note that dynamical simulations (AMD, BOB) performed for Xe+Sn 
central collisions at similar incident energies [t^, [^H predict average values 
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Fig. 30. Average excitation energy of primary fragments as a function of the atomic 
number for 50 MeV/nucleon Xe+Sn central collisions. Full points refer to exper- 
imental evaluations using fragment-particle correlations and open points to AMD 
simulations. From [72 



close to 3 MeV/nucIeon for fragment excitation at freeze-out. However a de- 
tailed comparison between data and AMD simulations, displayed in fig. [301 
shows different evolutions with fragment charge; the increase of the fragment 
excitation energy per nucleon with its charge is a general trend of AMD re- 
sults in the Fermi energy region (see also [8l|)- Going more into the details 
of correlation results, one learns that the corresponding secondary evaporated 
light charged particles represent less than 40% of all produced charged par- 
ticles and decreases down to 20% for 45-50 MeV/nucleon incident energies. 
Finally experimental data are better reproduced with secondary deexcitation 
simulations assuming that the fragment N/Z ratio at freeze-out is the same 
as that for the combined projectile-target system. 



1.3 Break-up densities and freeze- out volumes 



From an experimental point of view the location of nuclear multifragmenta- 
tion data in the phase diagram requires accurate independent measurements 
of temperature and density at the break-up stage. While the problem of tem- 
perature has been solved with acceptable accuracy using He-Li isotope ratios 



up to 5-6 MeV [15[, no accurate enough method is available to determine 
the spatial extension of mononuclear systems which undergo multifragmenta- 
tion. Therefore estimates of break-up densities have been obtained using var- 
ious approaches. For instance, break-up densities for spectator fragmentation 
in Au+Au collisions at 1000 MeV/nucleon were estimated by using selected 
particle-particle correlations (pa rticl es from secondary decays are excluded by 
imposmg an energy threshold) [iQs) . Assuming a zero lifetime, the volumes 
of spectator sources were extracted and densities calculated by dividing the 
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number of spectator constituents by the source volume. The estimated aver- 
age values slowly decrease from about 0.3 to 0.2 po when excitation energies 
of spectators move from 4 to 10 MeV/nucleon. Caloric curves were analyzed 
within the framework of an expa nding Fermi gas hypothesis to extract esti- 
mates of break-up densities 19d |. In this approach the observed flattening of 
the caloric curves reflects an increasing expansion with increasing excitation 
energy (the product Txp~^/^ is a constant within a isoentropic hypothesis). 
For nuclei of medium to heavy mass, the derived density values vary from 0.7- 
0.6 to 0.4 Po when excitation energies increase from 4 to 6-8 MeV/nucleon. 
Starting from the deflnition of freeze-out given at the beginning of that section, 
one can compare break-up density /volume to freeze-out volume obtained from 
dynamical simulations as BOB. Taking the case of Xe+Sn central collisions 
at 32 MeV/nucleon (source excitation energy close to 6 MeV/nucleon), break- 
up density coming from simulations is 0.41po or 2.4Vo and the corresponding 
average freeze-out volume VpQ^^ is 2.8Vo with the deflnition of SMM (which 
corresponds to a sphere including all the centres of primary fragments) and 



6.IV0 for Vp^Q a sphere including all fragments 78|, [l97|. As we see VpQ ^ 



' F.O. 

well compares, at least for that system size, to the break-up volume. 
Average freeze-out volumes have been extracted flrst from comparisons with 
statistical models like SMM as shown for some examples in subsection 5.2. 
In general, experimental results can be well reproduced over a large range 
of excitation energy (from 2 to 10, 12 MeV/nucleon), by keeping the same 
VpQ ^^ typically 3 or 6V0. Recently a flrst attempt of estimating the freeze-out 
volume in a fully consistent way was done, by employing a simulation built 
event by event from all the available asymptotic experimental information 
(charged particle spectra, av erage and standard deviation of fragment veloc- 
ity spectra and calorimetry) jl98 |. Dressed excited fragments, which statisti- 
cally deexcite, and particles at freeze-out are described by spheres at normal 
density. Four free parameters are used to recover the data: the percentage of 
measured particles which were evaporated from primary fragments, the col- 
lective radial energy, a minimum distance between the surfac es of products 
at freeze-out and a limiting temperature for fragments (see jl98 | for more 



details). For Xe-f-Sn data (32 MeV/nucleon, central collisions) an average vol- 
ume 

V^P^;=7.6±2.0H 

was derived. However the deduced average excitation 
energy for primary fragments is close to 4 MeV/nucleon, signiflcantly higher 
that the value of 2.2 MeV/nucleon evaluated from fragment-particle correla- 
tions (see previous subsection). Note that the use of the widths of fragment 
velocity spectra in the comparison between data and simulation has shown 
that the introduction of a limiting temperature for fragments seems manda- 
tory. The strong sensitivity of the freeze-out volume to limiting temperature 
values was also demonstrated in the microcanonical multifragmentation model 
MMM as far as average properties of fragments are used to put constraints on 



volumes 199 



In conclusion we can say that working hypotheses and approximations are 
used to give semi-quantitative information on average break-up densities or 
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freeze-out volumes. Values remaining constant around 0.2-0.3po or slowly de- 
creasing from 0.6 to 0.3po over the excitation energy range 3-12 MeV/nucleon 
are derived. The restoration of freeze-out stage using all the asymptotic ex- 
perimental information appears as promising. 



8 Finite systems and phase transitions 



Phase transitions are universal properties of interacting matter and tradi- 
tionally they have been studied in the thermodynamical limit of infinite sys- 
tems. A phase transition occurs when a phase becomes unstable in given ther- 
modynamical conditions described with intensive variables like temperature, 
pressure .... However in physical situations, as the one encountered in the 
present studies concerning isolated finite systems like nuclei, the concept of 
thermodynamical limit can not apply. Extensive variables like energy and en- 
tropy are no more additive due to the important role played by the surfaces 
of particles and fragments which are produced. The entropy of the surfaces 
which separate the coexisting phases does not scale with the size of the sys- 
tem. The entropy per particle at equilibrium s=S(E)/N shows a convexity 
with a depth proportional to A^~^/'^ which is suppressed at the thermody- 
namical limit. Consequently phase transitions should be reconsidered from 
a more general point of view. An important theoretical effort started a few 
years ago to propose and discuss s ignatures of phase tran s itions in smal l sys- 
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systems have been defined in reference |205l | as systems where the linear di- 
mension is of the same order of magnitude as the characteristic range of the 
interaction. This is for example the case for atomic clusters but also for astro- 
physical systems, due to the long range of gravitation. Such systems are also 
classified as non-extensive systems. 

This section is divided into three parts. One is devoted to the possible dy- 
namics of phase transition for hot unstable nuclei. In a second part we shall 
discuss the thermostatistics involved and the associated relevant signatures. 
Finally critical behaviours early observed in multifragmentation data will be 
discussed in the context of finite systems. 



8. 1 Dynamics and spinodal instability 



We have noticed in previous sections experimental evidence for a radial ex- 
tra energy boost (radial expansion energy) associated to multifragmentation 
products. It can be attributed either to a dominant compression-expansion 
phase in central nucleus-nucleus collisions or to thermal pressure for more 
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gentle collisions: hadron-nucleus or semi-peripheral nucleus-nucleus collisions. 
The system might then reach densities and temperatures that correspond to 
spinodal instability and clusterization would ensue as the system seeks to sep- 
arate into the corresponding coexisting liquid and gas phases. 
In the last fifteen years major theoretical progresses have been realized to un- 
derstand and learn about spinodal decomposition in the nuclear context and 
a review can be found in reference 22d |. We shall first briefiy discuss what 



are the specificities of spinodal decomposition as far as infinite nuclear matter 
is concerned. Associated to negative compressibility the mechanically unsta- 




Fig. 31. Nuclear matter dispersion re- 
lation at 3 MeV temperature for two 
different densi ties; po is the normal 
density, (from 
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Fig. 32. Growth rates of the most unsta- 
ble modes for a spherical source with 200 
nucleons function of the multipolarity 
L and for two different central densities, 
(from H). 

ble spinodal region can be inv estigated by studying the propagation of small 
density fiuctuations 228 . 23Cll |. By analogy with optics, the nuclear dispersion 
relation can be calculated for different conditions of temperature and den- 
sity by introducing the Boltzmann-Langevin equation (see subsection I5.3.2p . 
Within the linear response theory framework, if one expands the solution of 
the Boltzmann-Langevin equation as / = /o + Sf, where /o is a solution of 
the Boltzmann equation and 6f the fiuctuating part, one finds the equation 
of motion dSf /dt = —iA46f + 6I[fo] at the leading order in Sf. The extended 
RPA matrix Ai represents the combined action of the effective field and of 
the average collision term. In the spinodal region the eigenvalues of the matrix 
Ai become imaginary. Consequently the fiuctuations associated with a given 
eigenmode agitated by the source term 51 are exponentially amplified or sup- 
pressed, depending upon the sign of the imaginary part of the frequency. In 
the case of infinite nuclear matter the eigenmodes of the linearized dynamics 
are plane waves, characterized by a wave number k and an imaginary eigenfre- 
quency, which is the inverse of the instability growth time. Figure [31] presents 
an example of nuclear dispersion relation at 3 MeV temperature for two dif- 
ferent densities po/2 and po/3. Imaginary RPA frequencies are reported as a 
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function of the wawenumber k of the considered perturbation. This disper- 
sion relation exhibits a strong maximum at a given wave number followed 
by a cut-off at large k values. This cut-off reflects the fact that fluctuations 
with wavelength smaller than the range of the force can not be amplified. The 
most unstable modes correspond to wavelengths lying around A ~ 10 fm and 
the associated characteristic times are almost identical, around 30-50 im /c, 



depending on density (po/2-po/8) and temperature (0-9 MeV) }228l . |231| . A 



direct consequence of the dispersion relation is the production of "primitive" 
fragments with size A/2 ^ 5 fm which correspond to Z 10. However this 
simple and rather academic picture is expected to be largely blurred by several 
effects. We do not have a single unstable mode and consequently the beating 
of different modes occurs. Coalescence effects due to the residual interaction 



between fragments before the complete disassembly are also expected |228 
and finally primary fragments deexcite by secondary decay. 
Does the signal discussed for nuclear matter survive (in final fragment par- 
titions experimentally measured) if we consider the case of a hot expanding 
nucleus which undergoes multifragmentation? First of all, the system pro- 
duced by the collision has to stay long enough in the spinodal region (^3 
characteristic time: 100-150 fm/c) to allow an important amplification of the 
initial fluctuations. Secondly, the presence of a surface introduces an explicit 
breaking of the translational symmetry. Figure [32] shows the growth rates of 
the most unstable modes for a spherical sou rce o f A=200 with a Fermi shape 
profile and for two different central densities 229[ |. The growth rates are nearly 
the same for d ifferent multipolarities L up to a maximum multipolarity L^ax 
(see also 232| ). This result indicates that the unstable finite system breaks 
into different channels. Depending on multipolarity L, equal-sized "primitive" 
fragments are expected to be produced with sizes in the range AF/2-AF/Lmax] 
Ap being the part of the system leading to fragments during the spinodal de- 
composition. One can also note that the Coulomb potential has a very small 
effect on the growth rates of unstable collective modes except close to the 
border of t he s pinodal zone where it stabilizes very long wave-length unsta- 
ble modes 233|. On the other hand, for a finite system. Coulomb interaction 
reduces the freeze-out time and enhances the chance to keep a memory of 
the dynamical instabilities; a similar comment can be made as far as collec- 
tive expansion is concerned. Both effects push away the "primitive" fragments 
and reduce the time of their mutual interaction. So finally, even if expected 
strongly reduced, the presence of partitions with nearly equal-sized fragments 
is a good candidate to sign the role of spinodal instability in multifragmenta- 
tion. 



Following early studies related to nearly equal-sized fragment p artiti ons [234 



ten years ago a method called higher order charge correlations |235| was pro- 



posed to enlighten any extra production of events with specific fragment par- 
titions. The high sensitivity of the method makes it particularly appropriate 
to look for small numbers of events, as those expected to have kept a mem- 
ory of spinodal decomposition properties. Thus the charge correlation method 
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allows to examine model independent signatures that would indicate a pre- 
ferred decay into a number of equal-sized fragments in events from experi- 
mental data or from simulations. All fragments of one event with fragment 
multiplicity Mf = M = J2z''^z, where nz is the number of fragments with 
charge Z in the partition, are taken into account. By means of the normalized 
first, (Z), and second order, cr|, moments of the fragment charge distribu- 
tion in the event, one may define the higher order charge correlation func- 
tion as: 1 + R{az, {Z)) = {Y{az, {Z))/Y'{az, {Z)))\j^j. Here, the numerator 
Y{az, (Z)) is the yield of events with given (Z) and az values. The denomina- 
tor Y'^cTz, {Z)), which represents the uncorrelated yield of pseudo-e vents, can 
be built in different ways; details and discussions can be found in |72l . l236l . 12371 . 
The exact identification of fragments up to at least Z=20 is mandatory to use 
that method. It was applied on different data sets corresponding to com- 
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Fig. 33. Extra production of events with nearly equal-sized fragments (a: az < 1 and 
b: az < 3) as a function of thermal excitation energy (full points); the incident and 
radial energy scales are also indicated, e^^ and e^ad are deduced from comparisons 
with SMM. The open point refers to the result from BOB simulations; the average 
thermal excitation energy is used. Vertical bars correspond to statistical errors and 
horizontal bars refer to estimated uncertainties on the backtraced quantity, ej'^. 
(from[72]). 

pact single sources p rodu ced in central collisions in the incide nt en ergy range 
30-50 MeV/nucleon [isj. For the ^^^Xe+ """^Sn system 0, H a rise and 



fall of the percentage of extra production of events with nearly equal-sized 
fragments was observed. The average-charge domains contributing to the cor- 
relation peaks are 12-21 and come mainly, as theoretical ly ex pected, from 



fragment multiplicity 3 (dominance of the octupole mode) |239l | . Figure [33] il 



lustrates the results for two different limits on a^; the lower limit corresponds 
to the spread estimated as coming from secondary fragment evaporation and 
the higher limit to the observed spread for the correlation function. Results 
obtained with BOB simulations (see section [5]) are also reported. For the con- 
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sidered system, incident energies around 35-40 MeV/nucleon appear as the 
most favourable to induce spinodal decomposition; it corresponds to about 
5.5-6 MeV/nucleon thermal excitation energy associated to a very gentle ex- 
pansion energy around 0.5-1 MeV/nucleon. The qualitative explanation for 
those numbers can be well understood in terms of a necessary compromise 
between two time scales. On on e hand the fu sed systems have to stay in the 



spinodal region ^ 100-150 fm/c |228l . l23ll . |240| . to allow an important amplifi- 



cation of the initial fluctuations and thus permit spinodal decomposition; this 
requires a not too high incident energy, high enough however for multifragmen- 
tation to occur. On the other hand, for a finite system. Coulomb interaction 
and collective expansion push the "primitive" fragments apart and reduce the 
time of their mutual interaction, which is efficient to keep a memory of "prim- 
itive" size properties. 

At present confidence levels around 3-4 cr, observed for charge correlation 
peaks, prevent any definitive conclusion. To firmly assess or not the validity 
of this fossil signal new experiments with higher statistics are required. 
If this fossil signal is definitively confirmed the following comments have to 
be done. As we have shown in fig. [211 good average statistical descriptions 
of data using for example SMM are also obtained, which could demonstrate 
that the dynamics involved for finite systems is sufficiently chaotic to finally 
explore enough of the phase space in order to describe fragment production 
through a statistical approach. It is tempting to associate a part of that chaos 
with the coalescence stage which can occur during fragment formation. In- 
deed, as shown previously (sections [Hand [7]), the peculiar size and position of 
the largest fragment could be well understood as resulting from mode beat- 
ing (size), large- wavelength instabilities (size) and the late-stage coalescence 
related to the involved bulk nuclear density (size and position). 



8.2 Spinodal instability and Thermo statistics 



Spinodal instability is intimately related to the occurrence of a first order 
phase transition as signaled by a convex anomaly in the entropy function 
S{X). For pedagogical purpose we consider in the following only one single 
extensive variable X like energy, volume or number of particles. Entropy can 
be gained by separating the uniform system into two phases at equilibrium: 
it is the spinodal decomposition well known at the thermodynamical limit 
for binary solutions and binary alloys. Figure [31] illustrates the relationship 
between entropy convexity, spinodal instability and phase coexistence. Two 
signatures of a phase transition for hot nuclei can be added to the previous one 
(observation of nearly equal-sized fragments), just starting from the fact that 
surfaces are no more negligible for finite systems. The entropy of the combined 
system constituted by the two coexisting phases is not a sum of the individual 
subsystem entropies. The uniform unstable finite system gains entropy by re- 



50 




Fig. 34. Phase coexistence and spinodal instability. The entropy function S{X) is 
convex between the inflection points Xmin and X^ax where the intensive conjugate 
variable X{X)=6{S) /5{X) has a local minimum or maximum, respectively. The 
infinite system is mechanically unstable and entropy is gained by separating the 
system into two phases. The entropy of the corresponding mixed system is additive 
and moves along the common tangent as X increases from Xi to X2; the intensive 
conjugate variable is determined by th e Maxwell construction requiring that the 
two hatched areas are equal, (from 22d |). 



organizing itself into a mixture of the two phases but the resulting equilibrium 
entropy function keeps a local convexity. The Maxwell construction is no more 
valid as illustrated in fig. [351 

The two main consequences of that local convexity are the following. Using a 
microcanonical sampling (fixed values of the extensive variable E which can 
be rather well estimated from experiments) the first derivative corresponds 
to the inverse microcanonical temperature and the bottom curve of fig. [3l] is 
the inverse of a caloric curve which presents a backbending. The energetic cost 
paid to produce surfaces is a decrease of the microcanonical temperature when 
the excitation energy increases and con sequently the appearance of a negative 
microcanonical heat capacity 20Cll . 202] as a specific signature of a first order 



phase transition for finite systems. Moreover if a signature of spinodal insta- 
bility is observed one must measure correlatively a negative microcanonical 
heat capacity related to the resulting equilibrium entropy function with local 
convexity. The reverse is not true if spinodal instabilities are not responsible 
for fragment formation. Experimental observations of correlated signals will 
be discussed in the next section. 

If one considers now a finite system in contact with a reservoir (canonical sam- 
pling), the value of X may fluctuate as the system explores the phase space; 
the associated distribution at equilibrium is P(X)~exp(5'(X)-AX). The dis- 
tribution of X acquires a bimodal character (see fig. [36]): that bimodahty is 
a second specific signature of a first order phase transition for finite systems. 
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Fig. 35. Isolated finite systems. The en- 
tropy function for a uniform unstable 
system (lower curve) has a local con- 
vexity region and the system gains en- 
tropy reorganizing itself into a mixture 
of the two phases. The resulting equi- 
librium entropy function (upper curve) 
will always lie below the common tan- 
gent (dashed line) (from [22Cl[ |). 



Fig. 36. Canonical ensemble of 
finite systems. The bimodal equi- 
librium distribution is given by 
P(X)~exp(5(X)-AX). The figure 
shows the case when the Lagrange 
multiplier A is equal to the s lope of the 
common tangent (from 2201 ] ). 



It gives an understanding of coexistence as a bimodality of the event distri- 
bution, each component being a phase. We recall that for infinite systems 
two delta function distributions are localized in Xi and X2 and correspond 
to spinodal decomposition. The biggest fragment detected in each multifrag- 
mentation event, because of its correlation with the density which is the nat- 
ural order parameter of the liquid-gas phase transition, is a good candidate 
as a potential order parameter. A priori that second specific signature can 
be directly observed from experiments and could appear as robust. A diffi- 
culty comes however from the absence of a true canonical sorting in the data. 
The statistical ensembles produced by selecting for example fused systems 
are neither canonical nor microcanonical and should be better described in 



terms of the Gaussian ensemble [2411 ]. which gives a continuous interpola- 



tion between canonical and microcanonical ensembles. Very recently a simple 
reweighting of the probabilities associated to each excitation energy bin for 
quasi-projectil e eve nts was proposed to allow the comparison with the canon- 



ical ensemble |242l | . Signals of bimodality will be discussed and correlated to 
other signals in the next section. Finally one can also underline a very impor- 
tant consequence of bimodality: for finite systems large fluctuations of "order 
parameters" are expected to be observed in the middle of the coexistence region. 
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8.2.1 Caloric curves, negative microcanonical heat capacity and abnormal 
energy fluctuations 



The observation of a plateau in nuclear caloric curves was experimentally pro- 



178 . However, 



posed as a direct signature of a first-order phase transition 
from a theoretical point of view, a plateau-like shape can not be an unam- 
biguous signature even if it is a strong indication of a physical change and 
if its observation can help to better define the energy domain of interest for 
phase transition. Caloric curves for restricted mass regions are presented in 
fig. [271 Indeed measured caloric curves can be misleading because different 
curve shapes can be generated depending on the path followed in the mi- 
crocanonical equation of state surface. As examples calculated caloric curves 
(microcanonical lattice gasmodel-216 particles) at a constant pressure or a 
constant average volume |137l | are displayed in the upper part of fig. [371 At 
constant pressure a backbending is clearly seen whereas at constant average 
volume a smooth behaviour is observed showing that the phase transition sig- 
nal can be hidden in the observation of caloric curves. In experiments one does 
not explore a caloric curve at constant pressure nor at constant volume, the 
system follows a path in the excitation energy freeze-out volume plane. For 
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Fig. 37. Thermodynamic quantities in the microcanonical ensemble for a transfor- 
mation at constant pressure (left part) and at constant volume (right part). Caloric 
curves are displayed in upper panels. Normalized kinetic energy fluctuations are 
compared to the canonical expectation (lines) in lower panels. In medium panels 
microcanonical heat capacities (sy mbo ls) are compared to the estimation (dashed 



lines) discussed in the text, (from [137l ]). 



overcoming that situation it was proposed to directly measure the heat capac- 
ity; if an experimental path is intersecting the spinodal region of the system 
phase diagram, it would be signaled by a negative value of the microcanonical 
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heat capacity. A method for measuring m icrocanon i cal h eat capacity using 
partial energy fluctuations was proposed 137 . 204 . 24^. Abnormal partial 



energy fluctuations (as compared to the canonical expectation) are indeed al- 
ways seen, independently of the path, if microcanonical negative heat capacity 
is present. The prescription is based on the fact that for a given total energy 
of a system, the average partial energy stored in a part of the system is a 
good microcanonical thermometer, while the associated fluctuations can be 
used to construct the heat capacity. From experiments the most simple de- 
composition of the total energy E* is in a kinetic part, E^^, and a potential 
part, Ep, (Coulomb energy + total mass excess). However these quantities 
have to be determined at freeze-out and consequently it is necessary to trace 
back this configuration on an event by event basis. As discussed in the previ- 
ous section the fragment properties entirely rely on the representation of the 
system at the freeze out stage as non interacting fragments. The true config- 
uration needs the knowledge of the freeze-out volume and of all the particles 
evaporated from primary hot fragments including the (undetected) neutrons. 
Consequently some working hypotheses are used, possibly constrained by spe- 
cific experimental results [244] . Then, the experimental correlation between 
the kinetic energy per nucleon E^/A and the total excitation energy per nu- 
cleon E*/A of the considered system can be obtained as well as the variance of 
the kinetic energy cr^. Note that Ek is calculated by subtracting the potential 
part Ep from the total energy E* and consequently kinetic energy fluctuations 
at freeze-out reflect the configurational energy fluctuations. An estimator of 
the microcanonical temperature of the system can be obtained by inverting 
the kinetic equation of state: 

M o 

< Ek >=< ai > T'+ < -(M - 1) > T 
i=i ^ 

The brackets <> indicate the average on events with the same £"*, is the 
level density parameter and M the multiplicity at freeze-out. An estimate of 
the total microcanonical heat capacity is extracted from the following equa- 
tions: 



The specific microcanonical heat capacity {C j A)raicro becomes negative if the 
kinetic/configurational energy fluctuations Ao\ overcome CkT"^- Figure \37\ 
(medium and lower panels) illustrates the results of such a procedure in the 
framework of a microcanonical lattice gas model. 

That procedure was applied by the MULTICS and INDRA collaborations on 
different sources produced in peripheral an d central heavy ion collisions i n th e 



incident energy range 30-80 MeV/nucleon [158|, l245|, l246|, l247|, 12481, l249|, |250 



Figure [38] summarizes the results obtained by the MULTICS collaboration: on 
the left side it is seen that normalized kinetic energy fluctuations overcome 
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the estimated canonical fluctuations, Ck', the right part of the figure illustrates 
the microcanonical negative heat capacities observed, the distances between 
the poles being associated with the latent heat. 

These results provide a direct evidence of a first order liquid-gas phase transi- 
tion. They have to be seen as semi-quantitative and correspond to a starting 
point for more precise measurements needed for the reconstruction of the 
phase diagram. We recall again that the present protocol of measurements is 
discussed in details in reference 244 . Exact microcanonical formulae, assum- 






1 


i 


V 


- • 


■ ^ *l 




f 




1 

1 

1 


^ 





2.5 5 7.5 

EVAo (A.MeV) 



Fig. 38. Left panel: normalized kinetic energy fluctuations and estimated canonical 
fluctuations (C^) for quasi-projectile (QP) events produced in Au-|-Au collisions 
at 35 MeV/nucleon (grey zones) and for fused systems produced in central Au-|-C 
(black dots), Au-|-Cu (squares, triangles) and Au-|-Au reactions before (open stars) 
and after subtraction of 1 MeV/nucleon radial flow (black stars). Right panel: cor- 
responding microcanonical heat capacities per nucleon. (from (245j]). 

ing that a classical treatment of th e mo tion of products emitted at freeze-out 
is appropriate, are also proposed in 214| to calculate heat capacity and second 
order derivative of the entropy versus energy. They depend on the total kinetic 
energy and number of emitted products which have to be estimated event by 
event at freeze out. However up to now this method was not used to derive 
information from data. 



8.2.2 Microcanonical heat capacity and open questions 

Before concluding on this part one can indicate that theoretical questions are 
still under debate and concern the topology of the system at freeze-out and 
the related order of the phase transition. If the system is still relatively dense 
at freeze-out, which seems improbable if we refer to the previous section, the 
fragment properties may be very different from the ones asymptotically mea- 
sured and the question arises whether the energetic informatio n m easured on 
ground state properties can be taken for the freeze out stage |251 |. Classical 



molecular dynamics calculations have shown that the ground state Q-value 
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is a very bad approximation of the interaction energy of fragments in dense 
systems. This is due to the deformation of fragments and to the interaction 
energy when fragment surfaces touch each other. As a consequence, compa- 
rable kinetic energy fluctuations are obtained in t he su bcritical and super- 



critical region of the Lennard- Jones phase diagram [223|. On the other hand 



calculations with a similar model, the Lattice Gas model, show that even in 
the supercritical region the correct fluctuation behaviour can be obtained if 
both the total energy and the interaction energy are consistently estimated 
with the same approximate algorithm as it is done in the experimental data 



analysis [22J]. Concerning now the order of the transition, recent calculations 
performed within a mean-field framework (liquid drop model at finite temper- 
ature) show that the phase transition is seen to be continuous (second order) 
an d th at a negative heat capacity for finite nuclei is not incompatible with 
it [227) . 



8. 3 Scaling and fluctuations for fragment sizes 



In the early eighties multifrag mentation was connected to a cri t ical phe- 
nomenon 0, H, H, liil ^5M. Percolation models |25^, H, H, [issj, the 



Fisher Droplet Model |221l, |226|. l259l. |2 60. 261j and more recently the theory 



of universal fluctuations 20 tI . 2251 . 2621 . i263i . |264 | have been employed to de- 
scribe fragment size distributions and fluctuations and to tentatively derive 
information on the critical region of the liquid-gas transition: in finite systems 
the critical point becomes a critical region. The determinations of a scaling 
function and of a consistent set of critical exponents i n different multifragmen- 
tation data also tend to support this hypothesis [97|, Il4ll . |265| . |266| . However 



the situation appeared as more complex. Indeed signals of a first order phase 



97, 245, 266, 267 



transition were also observed on the same set of data 
which opened a debate on the critical temperature extracted, on the order of 
the transition and on the consequences of Coulomb and finite size effects on 



critical behaviours 203, 268, 269, 270, 271 



8. 3. 1 The Fisher Droplet Model 

The observation of critical behaviours like power laws in the charge/mass 
distribution of multifragmenting systems has been interpreted as a n evidence 



of a second order phase transition. In the Fisher Droplet Model |259l | the vapor 
coexisting with a liquid in the mixed phase of a liquid-gas phase transition is 
schematized as an ideal gas of clusters, which appears as an approximation 
for n on id eal fluids. The model was recently applied to multifragmentation 
data 260| by considering all fragments but the largest as the gas phase, Z^ax 



being assimilated to the liquid part. The yield of a fragment of mass A reads: 
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dN/dA = ri{A) = qoA-^ exp{{AAfx{T) - co{T)eA'') /T). 

In this expression, r and a are universal critical exponents, A/i is the difference 
between the liquid and actual chemical potentials, Co(T)eA'' is the surface free 
energy of a droplet of size A, cq being the zero temperature surface energy 
coefficient; e = {Tc — T)/Tc is the control parameter and describes the distance 
of the actual to the critical temperature. At the critical point Afi = and 
surface energy vanishes: rjlA) follows a power law. Away from the critical point 
but along the coexistence line A/x = and the cluster distribution is given by: 
dN/dA = r]{A) = qoA-^ ex.p{{-co{T)eA'')/T). 

The temperature T is determined by assuming a degenerate Fermi gas. This 




Fig. 39. Fragment mass yield distribution scaled by the power law prefactor, the 
chemical potential and Coulomb (added to AA/j, for left and central panels) terms 
plotted against the inverse temperature scaled by Fisher's parametrization of the 
surface energy. Left and central panels: data correspond to projectile fr agme ntation 
observed in reverse kinematics reactions at 1.0 GeV per nucleon, from 2611 ]. Right 
panel: Data correspond to semi-peripheral (Ta+X) or central (Xe+Sn, Ni+Au) 
collisions. Bombarding energies in MeV/nucleon are reported on to p of each scaled 
distribution; they are shifted by one x unit for a better view, from 



kind of scaling was found in many multifr agna e ntat i on data, in hadron-nucleus 



as well as in nucleus- nucleus collisions |260|, l26ll . l266l . l267l |; the agreement 
between data and theory often holds over orders of magnitude, and the critical 
exponents which are deduced are in acceptable agreement with those expected 
for the hquid-gas universality class. As examples, fig. [39] shows the scaling 
properties observed for different systems studied by the EOS and INDRA 
collaborations. 



8.3.2 Universal fluctuations of an order parameter 

The recently developed theory of universal scaling laws of order-parameter 
fluctuations provides methods to select order parameters, characterize critical 
and off-critical behaviours and determine critical point without any equilib- 
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rium assumption 2071 . |262| | . In this framework, universal A scaling laws of one 
of the order parameters, m, should be observed: 
(m)'^P(m) = 0((m — (m))/(m)^) 

where (m) is the mean value of the distribution P{m). A=l/2 corresponds to 
small fluctuations, cr^ ~ (m), and thus to an ordered phase. Conversely A=l 
occurs for the largest fluctuations nature provides, cr^ ~ ("^)^, in a disordered 
phase. For models of cluster pr oduc tion having a critical behaviour there are 
two possible order parameters [26^]: the fragment multiplicity in a fragmen- 
tation process or the size of the largest fragment in an aggregation process 
(clusters are built up from smaller constituents). The method was apphed 
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Fig. 40. (a) Zmax distributions for central Xe+Sn collisions at 25-39 MeV/nucleon 
bombarding energies, scaled according to A scaling equation; the dashed curve is a 
best fit to scaled data using a Gaussian distribution, (b) As (a) but for bombarding 
energies 39-100 MeV/nucleon: the d ashed curve is a best fit to scaled data using 
the Gumbel distribution. From |264 1 . 



to central collision samples (symmetric systems with total masses ~73-400 at 
bombarding ener gies between 25 and 100 MeV/nucleon) by the INDRA collab- 



oration 



263 



264l |. The CP multiplicity fluctuations do not show any evolution 
over the whole data set. Conversely the relationship between the mean value 
and the fluctuation of the size of the largest fragment does change as a func- 
tion of the bombarding energy: A ~l/2 at low energy, and A ~1 for higher 
bombarding energies. The form of the Z^ax distributions also evolves with 
bombarding energy: it is nearly Gaussian in the A=l/2 regime and exhibits 
for A=l an asymmetric form with a near-exponential tail for large values 
of the scaling variable (see fig . HOll . This distribution is close to that of the 
modified Gumbel distribution the resemblance increasing with the total 

mass of the system studied and being nearly perfect for the Au-|-Au data. The 
Gumbel distribution is the equivalent of the Gaussian distribution in the case 
of extreme values: it is obtained for an observable which is an extremum of a 
large number of random, uncorrelated, microscopic variables. 
The dependence on bombarding energy and total system mass of the frontier 
between the two A-scalings is shown in fig. HH Within the developed theory, 
this behaviour indicates the transition from an ordered phase to a disordered 
phase in the critical region, the fragments being produced following some ag- 
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Fig. 41. Dependence on bombarding energy and total system mass of the two differ- 
ent regimes of A-scaling observed for the size of the largest fragment in each event , 



Zmax^ for central collisions of symmetric systems studied with INDRA. From [2641 ] . 



gregation scenario. Here again, as it was the case for the Fisher's scaling, the 
two A-scaling regimes suggest the occurrence of a second-order phase transi- 
tion. 



8.3.3 Finite size effects and pseudo-critical behaviours inside the coexistence 
region 

As shown in this section, analyses of given samples of events have found at 
the same time scaling laws and evidence for first-order phase transition sig- 
nals (negative microcanonical heat capacity and spinodal instabilities). At the 
thermodynamic limit, scalings are general properties of matter near the criti- 
cal point and thus typical of second order phase transition. For finite systems 
fluctuations correlate points at a distance which can be comparable to the 
linear size of the system and can have an effect similar to a diverging corre- 
lation length in an infinite system. Therefore the so called critical behaviours 
are expected to be observed non only in the critical zone, but also deeply in 
the coexistence region where large fluctuations have been observed. Within a 
lattice gas model in the canonical ensemble a scaling behaviour, with criti- 
cal parameters close to the ones expected for the hquid-gas un ivers ality class. 



was observed for finite systems inside the coexistence region |203l |. Moreover 
a good quality scaling of the cluster size distribution was obtained using the 
Fisher's formula for the scaling function at subcritical as well as supercritical 
densities with values for the critical exponents comp atible (within finite size 
effects) with the universality class of the model 268| . This indicates that the 
observation of scaling does not allow to infer the position of the critical zone 
as it is also compatible with a fragmentation inside the coexistence region of a 
first order phase transition. Concerning the order parameter fluctu ation s, sim- 



ulations have been performed in the framework of the Ising model |273[ | which 
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show that the distribution of the heaviest fragment approximately obeys the 
A=l scahng regime even at subcritical densities where no continuous transi- 
tion takes place. The observed behaviour was interpreted as a finite size effect 
that prevents the recognition of the order of a transition in a small system. 
More recently the distribution of the heaviest fragment was analyzed within 



the lattice gas model |225[ | and it was shown that the most important finite 
size effect comes from conservation laws, the distribution of the order param- 
eter being strongly deformed if a constraint is applied (mass conservation) to 
an observable that is closely correlated to the order parameter. Moreover the 
observation of the A=l scaling regime was indeed observed in the critical zone 
but was also confirmed at subcritical densities inside the coexistence region. 

Very recently a generalization of Fishe r's model to extend and quantify fi- 



nite size effects was also proposed [226| . The Fisher's formula along the co- 



existence line is modified to the case of a vapor in equilibrium with a finite 
liquid drop replacing the original formula by dN/dA = r]{A) = go[^(^d " 
A)/A]-^exp{{-co{T)e)/T[A'^ + {A, - Af - A^]). 

Such a generalization was tested within the canonical lattice gas model. Bet- 
ter results are obtained for d=2 (square lattice) than for d=3 (cubic lattice). 
Applications to experimental data are needed to evaluate that potential im- 
provement. To conclude that part one can say that, from both experiments 
and theory, many different pseudo-critical behaviours can be observed inside 
the coexistence region of a first-order phase transition as far as finite systems 
are concerned. 



9 Coherent experimental signals of phase transition 



As shown in the previous sections, signals of phase transitions were observed 
by many experimental groups for systems with different sizes and excitation 
energies. In most cases however only one signal was evidenced, which is not 
sufficient to convince that a phase transition does occur, because each of the 
proposed signals has some intrinsic weakness. The quality of the exclusive 
measurements nowadays permits to search for several different signatures of 
the transition on a given sample of events; indeed the concomitant observation 
of several signals strongly reinforces the hypothesis that a phase transition 
has occurred. Several groups are presently looking for different signals in their 
samples. For instance, and if we restrict to reactions below 100 MeV/nucleon, 
the group of Bologna studied QPs from Au+Au collisions at 35 MeV/nu cleon , 
and after having discovered the negative branch of the heat capacity [245], 



they showed also Fisher scaling, critical exponents [266| and more recently 
they examine the possible occurre nce o f a bimodal distribution (see below) of 
the charge of the largest fragment 27J]. The INDRA collaboration evidenced 
several signatures of phase transition for Au QPs, and for central collisions 
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between Ni and Ni, Ni and Au, Xe and Sn (see ref 2371 . |248| . l267l | . A rather 
good coherence of the different signatures was found in these cases. We chose to 
focus in this section on three sets of data for which the most extended analyses 
were reahzed. They cover different system sizes and colhsion centrahty. The 
findings will only be summarized, the reader is sent to the quoted references 
for details: 

1) light quasi-projectiles (A~36), from reactions between 47 MeV/nucleon 
^OAr and ^^Al, ^^Ti and ^^Ni, by Ma ei a/ H; 



2) heavy quasi-pr oiectiles f rom Au on Au reactions, by the INDRA/ALADIN 
collaboration [7^, Il58l . l275l | ; 

3) fused systems of mass around 200 formed in centra l coll isions between Xe 



and Sn, by the INDRA collaboration [15|, l72|, l263|, [264, 1275 



Besides the signals mentioned in the previous section. Ma et al proposed some 
other ones as possible signature s of phase transition or critical behaviour, 
and firstly the nuclear Zipf law 2761 . and references therein]. Zipf law, first 
introduced in linguistic to analyse the relative population of words in English 
texts, says that the frequency of a word is inversely proportional to its rank, 
n, in a frequency list (the integer rank n is equal to 1 for the most frequently 
used word). It was afterwards evidenced in many different scientific areas and 
thus it was suggested that the Zipf hierarchy is a fingerprint of criticality. In 
nuclear physics, investigations with a lattice gas model showed that cluster 
distributions follow a Zipf law, where the frequency is replaced by the average 
size (charge or mass) of the clusters: oc n~^, where n = 1,2 . . . for the 
largest, second largest . . . fragment of each partition; A was found equal to 1 
(the Zipf law is verified) at the phase transition temperature. Ma fitted the 
average values of the charges of the multifragmentation clusters (starting at 
Z=l) as a function of their rank with this type of law for different excitation 
energies: A is found to decrease with increasing energy. The Zipf law is verified 
(A = 1) for the energy bin 5-6 MeV/nucleon. Consistently with the formulation 
of the Zipf law, percolation models suggest that the ratio of the sizes of the 
second to the first largest fragment s, S„ = {Zmax2) / {Zmax) , takes the value 0.5 
around the phase separation point [277|. And indeed, in the sample studied by 
Ma, Sp exhibits essentially linear behaviours, with a change of slope around 
the energy of 5.2 MeV/nucleon, where it takes the value of 0.5. 



As mentioned in section [8], for finite systems, in the transition zone of a first 
order phase transition, the distribution of the order parameter presents a bi- 
modal behaviour in the canonical ensemble, due to an anomalous convexity 

-1-1 -J 1 • ,„ • * 1 -L T^i -i- '-L' T^himodality 

m the underlymg microcanomcal entropy. Ihe transition energy, ^, 
corresponds to the value for which the two bumps of the bimodal distribu- 
tion have the same height. In the case of multifragmentation, the size of the 
heaviest cluster, which is correlated with the total energy deposit and the 
system density /volume, appears as a natural potential order parameter (see 
also the universal scaling laws in sect. 18.3.21) . The bimodal character of the 
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distribution may however be hidden if the variable is constrained by a con- 
servation law. Bimodality was searched for by using different observables con- 
nected with the charge partition of the events; the chosen variable in cases 1) 
and 3), where a limited excitation energy range was covered, has the form: 
Pb = (Ez,>Zh,„ - Ez,<Zu^ Zi)/ Ez,>Zth, with Ziim and Zth equal to 4 and 
1 for case 1), and 13 and 3 for case 3). Pb can be viewed as the difference 
between the liquid and the gas densities. In case 2) the distribution st udied 



was firstly that of the asymmetry between the two largest fragments [158 

); in a second step an elaborated analysis was 
performed, which consisted in weighting the Z^ax distribution by a probabil- 
ity given by the correlated excitation energy distribution on the total range 



scanned |74 1242 



The signals studied for the three systems investigated, and the excitation 
energy at which they occur are listed in table [H these energies are those com- 
ing from the calorimetry performed, and may thus include not only thermal 
energy but also some collective component. Scalings, fluctuations of the con- 
figurational energy and of the largest fragment charge, sudden changes in the 
evolution of key variables, topological structure of the fragment partitions, 
charge correlations were examined. 

(1) For the light A ~36 system, in the Texas A&M experiment, both charged 
products and neutrons were measured and used for centrality selections. A 
completeness criterion (Zgp >12) was applied, and QP formed in central 
collisions were analysed. The QP excitation energy, derived by calorime- 
try, ranges between 2 and 10 MeV/nucleon with a mean value around 
4.5 MeV/nucleon. Such a distribution does clearly not reflect the ex- 
pected impact parameter dependence, meaning that at low energy all the 
possible QP partitions are not explored. In this context, all the variables 
characterizing a phase transition studied present a coherent behaviour, 
with a change in their evolution in the range 5-6 MeV/nucleon. Note how- 
ever that, while a maximum of the configurational energy fluctuations is 
observed, its value remains very low (0.3), far below the canonical ex- 
pectation of 3/2. The authors give thus no firm statement about the 
observation of a negative heat capacity. Their global conclusion is that 
they show a body of evidence suggesting a liquid-gas phase change in an 
equilibrated system at, or extremely close to, the critical point. 

(2) For Au quasi-projectiles a completeness criterion is also required and 
calorimetry measurements lead to an excitation energy distribution de- 
creasing from 2 to 11 MeV/nucleon. For these heavier QP, with a charge 
~68, a similar consistency between all the signals of phase transition 
emerges. The maxima of fluctuations (fig. H2|) and E^**'*'^^ are in the 
range 4-5 MeV/nucleon, well inside the domain where the heat capac- 
ity is negative. The changes in the fragment topology (inversion of vari- 
ation of {Zmax2}i E^*'""'^"'**^ ) occur near the second divergence of the 
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Table 1 

Excitation energies - in MeV/nucleon - at which different signals were observed for 
three systems (see text). The energy values were derived from calorimetry. Zg is the 
charge of the considered source. E^*^'*'^'' is the excitation energy corresponding to 
Tc in Fisher formula (see sect. [8] and see text for the other definitions of variables). 
The threshold for radial expansion energy is reported for information (see sect. H]) 



variable 


QP As ~36 


QP Zs ~68 


monosources Zg ~82 


-r^Fisher 
crit 


5-6 


4.2 


3.8-4.5 


/\ tjnpi 1 1 n cr 








max /isfj^/ 1 


4-0. 








5-6 










\9 5-5 51 


f— -6 51 


max rr r7 iry 

^max 1 ^ s 




4-5 


< 5 


max {Zmax2) 


6 


5 


4.5-6 


^zipf- A = 1 


5.6 


8.5 


7.5 


Sp = 0.5 


5.2 


8.5 and above 


3.2 - 6 


change slope Sp 


5.6 


4 




-r-\bimodality 


5.6 


[4.75:5.25] 


7.8 


spinodal 




[5:8] 


[5:9] 


threshold Erad 




~5. 


~4.5 



(3) 



heat capacity, located at ~5 MeV/nucleon, as does the maximum of the 
extra equal-sized fragment events, possible witnesses of a spinodal de- 
composition. The values of E^**'**^^ and the position of the energies at 
which the heat capacity diverges agree within 10% with those obtained 
from ano t her sample of events obtained with the Multics/Miniball de- 
vice 245l . 1266( 1 . Conversely, while the variable Sp changes slope in the 
same domain, it reaches a value of 0.5 - where it saturates - at a much 
higher energy (8.5 MeV/nucleon, see fig. HSj); consistently the Zipf law is 
verified at the same high value (only fragments with Z>3 were taken into 
account for the heavy systems). In that case, the observations plead in 
favor of the occurrence of a first order phase transition of the liquid-gas 
type. The transition energy is lower than that derived for the lighter QP; 
the Zipf law and related signals can not be taken as a sign of transition. 
What can be said for the monosources formed in central Xe+Sn collisions 
at 25, 32, 39, 45 and 50 MeV/nucleon, which are slightly heavier than 
the QP of case 2)? They were isolated by a detected-charge completeness 
criterion, and by selecting events with a flow angle greater than 60° (com- 
pact event shapes in velocity space). The average excitation energies de- 
termined by calorimetry lie between ~5 and ~10 MeV/nucleon. In agree- 
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Fig. 42. Normalized fluctuations clS ct func- 
tion of the excitation energy for Au QP 
(Au + Au 80 MeV/nucleon). Open circles 
correspond to fluctuations of the largest 
fragment normalized by the source size 
Zi/Zg and full squares refer to normal- 
ized c onflg urational energy fluctuations, 
(from [27^ 



Fig. 43. Evolution of the variable Sp 
versus the excitation energy for Au 
quasi-projectiles and Xe-|-Sn compact 
sources from central collisions, data are 
from the INDRA collaboration. 



ment with the previous case, the fluctuations of both the conflgurational 
energy and the largest fragment charge, and the average charge of the sec- 
ond largest fragment, present a maximum in the region 4-5 MeV/nucleon 
where lies also E^*^'*'^''. Only the second divergence of the heat capacity is 
seen in the explored energy range, at a slightly higher than for Au QP's. 
The change of A scaling and the maximum of extra equal-sized fragment 
events occur near the second divergence, while the transition energy in- 
dicated by the bimodality signal is slightly higher (however the answer 
in that case is somewhat dependent on the chosen value of Zum)- Note 
that the Z^ax distributions are never bimodal in monosources, due to 
the strong constraints on the energy. Finally the Zipf law is also verifled 
above the second divergence, while the variable Sp takes the value of 0.5 
inside the coexistence region, but at an excitation energy varying with 
the incident energy (flgHSj). Let us underline that monosources with close 
mass and charge were also formed with an asymmetric entrance channel, 
Ni-|-Au, and that negative heat capacities and signal of spinodal decom- 
position we re o bserved in an energy range consistent with that quoted 
for Xe+Sn 

The consistency of all signals for central collisions is however not as good 
as in the Au QP case. It must be noted that for QP's radial expansion 
energy appears just above the higher limit of the coexistence region. Con- 
versely the threshold of expansion lies inside the coexistence region for 
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the monosources and may have a different influence on the signals, par- 
ticularly those related to the fragment topology: indeed as explained in 
sect, m radial expansion prevents coalescence between the nascent frag- 
ments, favouring the observation of the spinodal signal. 

In conclusion of this section, it may be said that a large body of signals of 
phase transition has been obtained in the three cases examined here. The tran- 
sition energy is higher when the system is lighter. The radial energy seems to 
influence the energy at which some signals occur, leading to a better overall 
consistency for QP's than for monosources formed in violent collisions. Finally 
there is no indication that the Zipf law is a signal of phase transition or criti- 
cal behaviour for the two heavy systems studied by the INDRA collaboration. 
The excitation energy, Ezipf, at which the Zipf law is verified was however 
shown to vary with the minimal value of the ch arge of the considered frag- 
ments: raising this value shifts Ezipj downwards |275| . One may thus suspect 



a connection between the system size and the range of charge to consider in 
order to establish the Zipf law. While Ma et al. conclude that "the transition 
occurs close to the critical point" for their light system, for heavy systems 
the ensemble of results inclines towards the occurrence of a first order phase 
transition. 



10 Conclusions 



With the advent of powerful multidetectors dedicated to the study of mul- 
tifragmentation, which allow in particular the detection of more complete 
events, important improvements have been realized in the sorting of data 
and the construction of global observables. New correlation methods were 
developed (higher order charge correlations, fragment-particle) to investigate 
fragment formation mechanism and primary fragment excitation energy. Yield 
scaling as well as various fluctuation properties were also studied in relation 
with phase transition. Detailed and instructive comparisons with dynamical 
and statistical models have been reported. 

Fragment emission exhibits both dynamical and statistical aspects. For ex- 
ample, heavy-ion collisions and hadron-nucleus collisions were used to clearly 
show that fragment radial expansion energy is higher for central heavy-ion 
collisions for which a compression-expansion phase is predicted to occur in 
dynamical simulations; deduced radial expansion energies from experiments 
are around 4-5 MeV/nucleon at the maximum of fragment production. Frag- 
ment multiplicity appears as mainly governed by thermal energy whereas frag- 
ment partitions seem to be sensitive to radial expansion, which could be an 
indication of different trajectories in the phase diagram. Note that statisti- 
cal models predict the sensitivity of partitions to break-up density/volume. 
Fragment-particle correlation studies strongly suggest that primary fragments 
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are, on average, in thermal equilibrium. That result is supported by both 
the success of statistical models in describing rather well average multifrag- 
mentation properties and the information derived from dynamical simulations 
concerning the different time scales involved along the coUisions. A saturation 
of the average excitation energy of primary fragments around 3 MeV/nucleon 
is also deduced. Moreover a few recent examples show a rather coherent link 
between dynamical and statistical descriptions. 

What is the mechanism of fragment formation? It is still an open question. In 
stochastic mean field approaches the fragmentation process follows the spin- 
odal decomposition scenario whereas in molecular dynamics many-body corre- 
lations play a stronger role and pre-fragments appear at earlier times. New ex- 
periments with larger statistics are needed to eventually demonstrate through 
fragment charge correlations that spinodal instabilities could be responsible 
of fragmentation. On the theoretical side, although large progresses have been 
made, the development of quantal transport theory is still a challenge to for- 
mally describe the basic quantal nature of hot nuclear systems and of their 
fragmentation. 

The specific role of the largest fragment in each event was also clearly re- 
vealed. Firstly, at a given excitation energy, its charge distribution is inde- 
pendent of the total mass of the multifragmenting system when the latter 
exceeds A ~ 180-190; this behaviour is also observed in statistical models 
down to A~ 150. Then, its specific position close to the centre of mass of the 
multifragmenting systems produced in central heavy-ion collisions, and more 
randomly positioned for ^He-nucleus collisions, informs on event topology at 
freeze-out. Finally, its scaling properties and bimodal distributions can indi- 
cate that the largest fragment may be considered as a good candidate as order 
parameter of a phase transition in nuclei. 

Since the early eighties multifragmentation has been tentatively associated 
with the liquid-gas phase transition theoretically predicted for nuclear mat- 
ter at sub-saturation densities and temperature below 16-18 MeV. However, 
when dealing with the breaking of a nucleus, statistical mechanics of finite 
systems appeared as a key issue to progress and propose new first-order phase 
transition signatures related to thermodynamic anomalies like negative micro- 
canonical heat capacity and bimodalities. Those anomalies are generic features 
in non-extensive systems and have been reported in other fields such as atomic 
cluster physics and self-gravitating systems. 

A big effort to accumulate experimental indications of a phase transition has 
been presently made, including scaling and fiuctuations for fragment sizes. At 
the same time progresses have been obtained to experimentally better localize 
nuclear multifragmentation in the nuclear phase diagram thanks to improve- 
ments in calorimetry, thermometry and estimates of density at the break-up 
stage within working hypotheses and approximations. Today a rather coherent 
picture has been reached for few exhaustive studies and some semi-quantitative 
estimates of the transition region become available. However the order and the 
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nature of the transition is still subject to debate. For a large part this debate 
is related to finite-size effects, which remain an important challenge. 

How to progress further ? Experimentally new advances can be made by iden- 
tifying charge but also mass of fragments with the construction of a new gener- 
ation of 4tt multidetectors. With such new devices some hypotheses would be 
reduced or suppressed to evaluate freeze-out properties from intra-event cor- 
relations and it would become possible to disentangle radial expansion from 
Coulomb repulsion. On the theoretical side, by largely varying the proportions 
of neutrons and protons, a new phenomenon related to the phase transition for 
nuclei is predicted: the distillation, which makes the gas phase more asymmet- 
ric than the liquid phase (even more asymmetric for spinodal decomposition 
as compared to phase equilibrium). By measuring the isotopic composition 
of all the fragments it would become possible to identify the two phases and 
determine by a new way the dynamics of the transition. 
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